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EXECUTIVE  SUMMARY 


In  May  of  1979  the  Massachusetts  Department  of  Environmental  Quality 
Engineering  (DEQE)  discovered  that  two  wells  in  Woburn,  Massachusetts  were 
contaminated  with  toxic  chemicals.   Public  health  studies  have  indicated 
that  a  high  incidence  of  childhood  illness  has  occurred  in  the  general  area 
near  the  contaminated  wells.   A  study  of  Woburn's  water  distribution  system 
was  done  by  the  New  England  Interstate  Water  Pollution  Control  Commission  to 
define  the  area  affected  by  these  wells.   Since  the  definition  of  the  area 
affected  by  the  contaminated  wells  is  important  to  showing  the  connection 
between  the  wells  and  the  childhood  illness,  this  report  presents  a  second 
water  distribution  analysis  to  complement  and  evaluate  the  results  of  the 
first  water  distribution  study. 

The  method  used  to  develop  this  second  water  distribution  analysis 
begins  with  the  data  collection  and  calibration  of  a  computer  model  of  the 
present  (1983-1985)  distribution  system.   While  the  computer  program  used 
for  the  network  analysis  is  a  standard,  commercially  available  analysis 
procedure,  the  data  available  were  somewhat  unusual.   Good  data  for  the 
spatial  distribution  of  the  water  consumption  and  the  water  pipe  roughness 
were  not  available.   Therefore  data  from  fire  hydrant  tests,  done  for 
insurance  purposes,  was  used  to  calibrate  the  computer  model. 

The  computer  model  of  the  1983-1985  water  system  was  then  field  tested 
to  evaluate  the  accuracy  of  the  model.   Good  agreement  was  found  between  the 
measurements  and  the  computer  model. 

The  1983-1985  computer  model  was  then  modified  to  correspond  to  the 
historical  water  systems  of  the  period  when  the  contaminated  wells  were  in 

ix 


use.   Two  separate  historical  periods,  1964-1969  and  1970-1979  are  needed 
since  the  water  system  changed  during  the  pumping  history  of  the 
contaminated  wells.   Further,  the  on-off  status  of  one  of  the  non- 
contaminated  wells,  the  percentage  of  the  total  water  supply  that  came  from 
the  contaminated  wells,  and  the  daily  cycle  of  water  usage  all  affected  the 
spatial  distribution  of  the  contaminated  water.   Thus  the  distribution 
patterns  of  the  Woburn  water  system  were  found  to  be  quite  varied,  much  more 
varied  than  was  indicated  by  the  previous  study.   Figures  2H   to  35  present 
these  results. 

This  second  study  indicates  many  errors  in  the  first  study,  but  even 
with  a  large  number  of  significant  differences  between  the  two  analyses  the 
results  are  roughly  the  same.   Both  studies  indicate  a  southern  boundary  to 
the  contaminated  water  near  Montvale  Avenue  and  a  western  boundary  near  Main 
Street.   While  this  general  result  is  the  same  there  are  significant 
differences  in  the  results.   The  principal  difference  is  the  greater 
variation  in  the  location  of  the  boundaries  of  the  contaminated  water  that 
is  indicated  by  this  study. 


I.   Introduction 

In  May  of  1979  the  Massachusetts  Department  of  Environmental  Quality 
Engineering  (DEQE)  discovered  that  two  wells  in  Woburn,  Massachusetts  were 
contaminated  with  toxic  chemicals.   These  wells  were  shut  off.   The  rest  of 
the  city's  wells  were  not  contaminated  and  continue  in  operation.   The 
Metropolitan  District  Commission  (MDC)  now  supplies  the  water  required  to 
replace  the  water  from  the  two  wells. 

Studies  of  the  health  of  the  people  of  Woburn  showed  that  a  part  of 
the  city  had  a  higher  level  of  childhood  illness  than  would  normally  be 
expected.   The  affected  part  of  Woburn  roughly  coincided  with  the  area 
served  by  the  two  contaminated  wells.   The  period  when  the  excess  illness 
occurred  also  indicated  a  possible  connection  between  the  well  water  and  the 
public  health  problem. 

In  198-4  more  detailed  research  has  shown  "a  statistically  significant 
positive  association"  between  the  incidence  rate  of  various  childhood 
diseases  and  the  use  of  water  from  the  contaminated  wells.   The  water 
distribution  system  of  Woburn  was  analyzed  by  Helen  A.  Waldorf  and  Robert  K. 
Cleary  of  the  New  England  Interstate  Water  Pollution  Control  Commission. 
Based  on  their  1984  report  (6)  and  also  on  previous  estimates  of  the  area  of 
Woburn  supplied  with  water  from  the  contaminated  wells,  a  1984  journal 
article  (4)  by  S.  W.  Lagakos,  B.  J.  Wessen  and  M.  Zelen  has  indicated  that 
the  high  incidence  rates  of  childhood  illness  are  associated  with  the  area 
served  by  the  two  contaminated  wells. 

Since  the  definition  of  the  area  affected  by  the  contaminated  wells  is 
so  important,  a  second  water  distribution  analysis  was  proposed.   Thus  the 
present  study  evaluates  and  complements  the  study  of  Waldorf  and  Cleary. 


There  have  been  many  hydraulic  studies  of  the  water  system  of  Woburn 
during  the  historic  period  of  Weils  G  and  H.   In  1958  Whitman  and  Howard  (7) 
recommended  six  major  improvements  to  the  water  supply  system  including  the 
construction  of  the  Rag  Rock  reservoir.   This  study  was  quite  comprehensive, 
considering  the  user  demands,  ground  water  supply  and  the  status  of  the 
distribution  system  at  that  time.   In  1970  Coffin  and  Richardson  (1)  studied 
the  water  supply  needs  of  Woburn  and  other  towns  north  of  Boston  for  the 
MDC.   They  too  looked  at  user  demands  and  analyzed  the  flow  distribution, 
and  concluded  that  Woburn  needed  help  from  the  MDC  system.   They  recommended 
a  first  connection  to  the  MDC  network  at  Montvale  Avenue  and  Washington 
Streets  as  well  as  a  second,  later  connection  in  North  Woburn.   In  1970 
also,  Whitman  and  Howard  (8)  studied  the  water  needs  of  North  Woburn  and 
recommended  both  the  MDC  connections  and  the  installation  of  new  mains  in 
North  Woburn.   In  1980  the  Pipe  Line  Testing  Service  Company  performed  a 
water  leakage  study  (5)  to  determine  the  extent  of  any  leakage  in  the  water 
distribution  system.  Their  report  indicated  very  substantial  leakage 
totalling  1.7  million  gallons  per  day.   Although  there  is  some  doubt 
concerning  the  accuracy  of  this  amount,  it  was  clear  that  Woburn  had  a 
serious  leakage  problem.   Since  residential  users  are  billed  at  a  flat  rate 
rather  than  from  meter  readings,  there  is  only  a  poor  description  of  the 
areal  distribution  of  water  use  in  Woburn  at  the  present  time.   The  most 
recent  study,  that  of  Waldorf  and  Cleary  (6),  was  focused  on  the  present 
problem,  "Where  did  the  water  from  Wells  G  and  H  go?"  They  examined  the 
water  distribution  system  during  the  1964-1979  period  when  wells  G  and  H 
were  in  service  and  developed  two  separate  computer  models  for  the  system, 
one  for  the  1964-1969  period  and  one  for  1970-1979.   Based  on  analysis  of 


the  fraction  of  the  Woburn  supply  provided  by  wells  G  and  H  they  chose  four 
typical  months  for  each  period  and  calculated  the  flow  distributions  for 
each  month. 

As  was  mentioned  previously  the  goal  of  the  present  study  is  the  same 
as  that  of  Waldorf  and  Cleary,  to  determine  which  part  of  Woburn  was 
supplied  with  water  from  wells  G  and  H.  Thus  this  study  examines  the  data 
needed  to  describe  the  water  system,  develops  a  computer  model  for  the 
analysis  of  the  system,  validates  the  model  with  a  field  test  and  then  uses 
the  model  to  calculate  the  areas  supplied  with  water  from  the  contaminated 
wells. 
II.   Woburn  Water  System  Data 

Overall,  the  data  collected  is  more  than  sufficient  for  the  purpose  of 
the  study.   A  source  list  for  these  data  is  presented  in  Appendix  I. 

The  water  pipes  were  mapped  by  Thomas  Mernin,  Woburn  City  Engineer  and 
in  greater  detail  (1"  =  100')  by  Lockwood,  Kessler,  and  Bartlett.   Their 
maps  show  the  location,  length,  diameter,  material,  age,  and  interconnection 
of  the  Woburn  water  pipe  network. 

In  this  study  the  Woburn  water  pipe  network  was  simplified.   The 
principal  simplification  was  the  elimination  of  pipes  that  were  not  part  of 
pipe  loops.   The  user  demand  of  the  omitted,  "dead  end"  pipes  was  assumed  to 
occur  at  the  nearest  pipe  junction  of  the  remaining  "looped"  pipes.   The 
second  simplification  was  elimination  of  small  pipes  whose  carrying  capacity 
was  less  than  one  hundredth  of  larger  pipes  connecting  the  same  pipe 
junctions.   Typically  these  small  pipes  were  six  inches  in  diameter,  while 
their  larger  neighbors  were  16  inches  in  diameter.   Finally,  equivalent 
lengths  of  pipe  were  used  to  simplify  cases  when  two  different  diameters 


were  found  between  junctions  and  when  a  number  of  small  pipes  were  located 
in  parallel. 

Figures  1,2,  and  3  show  the  simplified  pipe  networks  of  Woburn  for  the 
1964-69,  1970-79,  and  1983-85  periods  respectively.   Table  1  presents  these 
data  also. 

There  were  two  problems  associated  with  the  pipe  data,  the 
interconnection  of  old  and  new  pipes  along  the  same  streets  and  the 
hydraulic  roughness  of  the  pipes.   The  first  problem  occurs  along  Cambridge 
Road,  Bedford  Road,  Winn  Street,  and  near  the  Horn  Pond  Reservoir. 
Interviews  with  the  water  personnel  of  the  City  of  Woburn  clarified  the 
interconnections. 

The  hydraulic  roughness  problem  is  more  difficult.   The  report  by 
Whitman  and  Howard  (7)  in  1958  describes  the  measured  roughness  of  eight 
typical  water  mains  and  that  description  is  presented  in  Table  2.  The 
average  C  value  of  this  sample  is  80.6  while  the  average  age  of  the  measured 
pipe  (in  1957)  was  34  years.   Since  it  is  generally  accepted  that  clean,  new 

pipe  starts  with  C  =  130  and  that  C  decreases  as  the  age  of  the  pipe 
increases,  Figure  4  presents  a  description  of  the  aging  process  that  fits 

the  Woburn  roughness  data.  This  figure  implies  that  C  =  40  for  those  pipes 
in  1985. 

Personal  examination  of  four  pieces  of  ten  inch  and  six  inch  cast  iron 
pipe  that  had  been  removed  from  various  Woburn  streets  showed  that  the  inner 
wall  roughness  of  the  pipe  was  nearly  1/8th  inch  for  all  samples.   A  single 
piece  of  asbestos-cement  pipe  had  a  roughness  of  1/32nd  inch.   Using  the 
Reynolds  Number  and  relative  roughness  of  these  pipes  to  estimate  the 
friction  factor,  f,  of  the  pipe  flows  permits  the  evaluation  of  the 
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Bedford 
Bedford  and  Hillside 
Kilby,  Hart,  and  Wyman 
Fowle  and  Garfield  to 

Green 
Green  from  Main  to  Garfield  1700 
Montvale  from  Main  to 

Auburn  1100 


1100 

10 

CI 

1910 

75 

33 

18 

27 

2100 

16 

CI 

1959 

26 

31 

22 

27 

1750 

20 

CI 

1959 

26 

35 

27 

28 

moo 

21 

CI 

1959 

26 

36 

28 

29 

3800 

16 

CI 

1959 

26 

37 

28 

26 

1500 

8 

CI 

1915/27 

70/58 

38 

19 

30 

1700 

10 

CI 

1929 

56 

39 

20 

30 

12 


CI 


1922 


63 


10 


21 


31 


Auburn  and  Garfield  to 

Green  3100 

Montvale  from  Auburn  to 

Bow 
Montvale  from  Bow  to  Wood 
Montvale  from  Wood  to 

Green 
Green  from  Garfield  to 

Montvale 
Montvale  from  Green  to 

Central 
Montvale  from  Central  to 

Washington 
Montvale  from  Washington 

to  MDC 
Central  and  D  Streets 
Washington  from  Montvale 

to  D  2500 


CI 


1927 


58 


CI 


1920 


65 


11 


30 


50 


36 


31 


1900 

12 

CI 

1922 

63 

12 

31 

32 

2200 

12 

CI 

1922 

63 

13 

33 

32 

1000 

12 

CI 

1910 

75 

11 

31 

33 

1500 

10 

CI 

1929 

56 

15 

31 

30 

1500 

12 

CI 

1910 

75 

16 

35 

31 

600 

12 

CI 

1910 

75 

17 

36 

35 

100 

12 

CI 

1975 

10 

18 

37 

36 

3300 

6 

CI 

1911 

71   • 

19 

35 

38 

38 


Washington  from  Montvale 

to  Pine 
Central  from  Montvale  to 

Pine 
Salem  from  Main  to  Bow 
Bow 

Salem  from  Bow  to  Wood 
Wood 

Salem  from  Wood  to  Pine 
Pine 
Washington  from  Pine  to 

Cedar 
Salem  from  Pine  to  Wells 

G  &  H 

Salem  and  Cedar  to 

Washington 
Wells  G  &  H  Connector 
Clinton 
Beach 
Eaton 


1800 


CI 


1932 


1300 

12 

CI 

1975 

3100 

10 

CI 

1920 

2000 

6 

CI 

1921 

1500 

10 

CI 

1920 

2600 

6 

CI 

1920 

1200 

10 

CI 

1920 

3500 

10 

AC 

1930 

3100 

6 

CI 

1930 

1700 

12 

AC 

1965 

1100 

12 

AC 

1965 

750 

12 

AC 

1965 

800 

6 

CI 

1928 

3000 

6 

CI 

1928 

1100 

16 

CI 

1960 

53 

10 
65 
61 
65 
65 
65 
55 

55 

20 


20 

20 
57 
57 
25 


51 


36 


39 


52 

35 

39 

53 

23 

10 

51 

32 

10 

55 

10 

11 

56 

33 

11 

57 

12 

11 

58 

39 

12 

59 

39 

11 

60 

12 

13 

61 

13 

11 

62 

15 

13 

63 

25 

16 

61 

10 

16 

65 

26 

17 

Table    1    (continued) 


Mishawum  from  Clinton 

to  Eaton 
Mishawum  from  Eaton  to 

Forest  Pk. 
Mishawum  from  Forest  Pk.  to 

Olympia 
Olympia  from  Mishawum  to 

Wi ldwood 
Wildwood 

Olympia  from  Wildwood  to 

Washington 
Washington  from  Cedar  to 

Olympia 
Forest  Pk.  from  Mishawum  to 

Alfred 
Forest  Pk.  from  Alfred  to 

School 
School  from  Forest  Pk.  to 

Mishawum 
Mishawum  from  Olympia  to 

School 
New  Boston  from  School  to 

New  Industrial 
New  Industrial 
Ryan 
Mishawum  from  Ryan  to 

Commerce 


600 

6 

CI 

1957 

28 

66 

46 

47 

600 

16 

CI 

1960 

25 

67 

47 

48 

2500 

16 

CI 

1960 

25 

68 

48 

49 

1000 

16 

CI 

1960 

25 

69 

49 

50 

5000 

12 

CI 

1970 

15 

70 

41 

50 

3800 

16 

CI 

1960 

25 

71 

50 

51 

2000 

12 

AC 

1965 

20 

72 

44 

51 

2800 

16 

AC 

1970 

15 

73 

48 

52 

1100 

16 

AC 

1970 

15 

71 

52 

53 

1600 

6 

CI 

1932 

53 

75 

53 

54 

1000 

6 

CI 

1899 

86 

76 

49 

54 

1100 

16 

AC 

1970 

15 

77 

53 

55 

1800 

12 

AC 

1971 

14 

78 

55 

56 

800 

12 

AC 

1971 

11 

79 

54 

56 

2500 


16 


CI 


1978 


80 


56 


57 


Mishawum  from  Commerce  to 

Washington 
Washington  from  Olympia  to 

Mishawum 
Commerce  Way 
Commonwealth  Avenue 
Washington  from  Mishawum 

to  Commonwealth 
New  Boston  from  New 

Industrial  to  Merrimac 
Main  from  Wyman  to  Alfred 
Alfred 

Elm  from  Main  to  West 
Elm  from  West  to  Main 

Main  from  Elm  to  Nichols 
School  from  Main  to 

Merrimac 
School  from  Merrimac  to 

New  Boston 
Merrimac 
Nichols 
Webster 
Pearl  and  We3t 


1100 


3300 


24 


CI 


CI 


1973 


12 


1920/38    65/47 


81 


85 


58 


58 


57 


1600 

16 

CI 

1973 

12 

82 

51 

58 

3200 

20 

CI 

1973 

12 

83 

57 

59 

3700 

10 

AC 

1977 

8 

84 

59 

60 

60 


4100 

16 

AC 

1970 

15 

86 

55 

61 

4100 

16 

CI 

1917 

68 

87 

26 

62 

2900 

6 

CI 

1931 

54 

88 

52 

62 

1300 

12 

CI 

1920 

65 

89 

62 

63 

1800 

12 

CI 

1920 

65 

90 

63 

64 

500 

12 

CI 

1927 

58 

91 

64 

65 

700 

6 

CI 

1924 

61 

92 

64 

66 

3700 

6 

CI 

1932 

53 

93 

53 

66 

4900 

6 

CI 

1924 

61 

94 

66 

61 

1200 

12 

CI 

1927 

58 

95 

65 

67 

1500 

12 

CI 

1927 

58 

96 

67 

68 

3000 

6 

CI 

1916 

69 

97 

63 

68 

10 


*CI   =  Cast   Iron 
AC  =   Asbestos   Cement 


Table    1    (continued) 


Pool  and  Winter 

2000 

6 

CI 

1925 

60 

98 

67 

69 

Mounta  in 

2200 

6 

CI 

1923 

62 

99 

70 

69 

Main  from  Nichols  to 

Mountain 

2000 

6 

CI 

1913 

72 

100 

65 

70 

11 


Table  2.   Pipe  Roughness  in  1958(7)  -  Cast  Iron 

Street  Diameter( in. ) 

Cambridge  6 

Merrimac  6 

Green  10 

Salem  10 

Bedford  10 

Montvale  12 

South  Bedford  12 

Arlington  16 


Date  Built 

"C" 

1921 

100 

1924 

85 

1921 

70 

1920 

65 

1937 

70 

1922 

75 

1919 

90 

1920 

90 

12 


140- 


120- 


100- 


U 


Whitman  &  Howard  (7) 
1958 


1920 


1940 


1960 


1980 


Date 


Figure  4 

Pipe  Roughness 
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corresponding   Hazen-Williams   C  values.      All   five  samples  were   assigned   C 
values  between  90   and   100  by   this  procedure.      Since   the  Woburn   water  workers 
said   these  five   were   typical  of   the  many   other  water  pipes   they   dug   out  of 
the   streets,    this   evidence  places  C  near   100   and  contradicts   the  previous 
data. 

Rather   than  'try   another   C  measurement  program   in   the  winter  of  1985, 
some   fire   hydrant   test   data  of  the  ISO  Commercial   Risk   Services   Company   will 
be  used   to   "adjust"   the  C  values   of   the  Woburn   pipes.      This   adjustment 
process   is   explained  later   in   Section   III  of   the  report. 

Only   a   few  key   valves    in   the   pipe  network   deserve   comment.      There    is   a 
set  of  valves  just   west  of   Cambridge   Street   that   isolate   the   Shaker  Glen 
high  pressure  sub-system  from  the  rest  of  the  Woburn  pipe  network.      The  only 
connection  between  the  pipes  west  of  Cambridge  Street  and  those  along  and 
east  of  it   is   through  the  booster  pump  at  Russell  and  Cambridge   Streets. 
Also,    there   is  a  closed  valve  on  Burlington   Street  near   Lexington   Street. 

The  water  billing   system  of  Woburn  provides   the  best   available 
information  on  the  spatial   distribution  of   the  city's   water   consumption. 
Water  meters,    read   twice  each  year,    record   the  amounts   supplied   to   large 
private  users.      In   1984   these  meters   accounted   for  21    percent  of   the   water 
supplied.      Large  municipal  users,    like   the  high  school,    are  metered   but   the 
meters   are  no   longer   read  since   the  city   would  be  billing  itself.      The  water 
department   and   water   service  are   not   an   isolated   financial  unit   of  the  city. 
Further,    not  all   the  residences  of  Woburn  are  metered.      Rather,    all 
residences   are   charged  a  flat-rate   annual   fee   for   water.      With   this 
residential  billing   system  there   is   little  reason   to  read   the  meters  of 
those   houses   that   have   them,    so  residential  meter   reading   has   been 
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discontinued.   However,  even  unmetered  water  bills  are  useful  since  they 
provide  the  addresses  of  residential  water  users. 

Coffin  and  Richardson  (1 )  reported  that  residential  and  commercial 
consumption  was  85  gallons  per  capita  per  day  during  1969  and  that  this  user 
demand  was  gradually  increasing  with  time  to  100  gallons  per  day  now  (1985). 
The  census  data  of  1960  shows  an  average  of  3.7  people  per  household  and  so 
a  "good"  estimate  of  370  gallons  per  day  per  residential  water  meter  is 
based  on  this  data. 

A  street  by  street  distribution  of  water  demand  was  prepared  from  the 
water  bills  of  1984.   This  distribution  was  then  simplified  by  subdividing 
Woburn  into  a  number  of  user  demand  areas  corresponding  to  the  water  supply 
pipe  junctions.   Each  user  street  was  assigned  to  is  nearest  junction  along 
the  path  of  the  water.   Figure  5  shows  a  map  of  the  user  demand  areas  and 
the  corresponding  pipe  junctions  or  nodes.   Table  3  indicates  the  boundary 
points,  along  the  pipes,  for  each  user  area.   Appendix  II  gives  the  street 
list  for  each  user  area.   The  residential  demand  of  each  user  demand  area 
was  added  to  the  metered  demand  of  that  area  and  the  resulting  sums  are 
given  in  Table  4.   The  combined  industrial  (1.29  MGD)  and  residential  (4.47 
MGD)  demand  accounts  for  95  percent  of  the  water  supplied  (6.06  MGD)  in 
1984. 

While  the  industrial  and  residential  users  were  easily  located  through 
the  water  billing  system,  the  spatial  distribution  of  the  water  leakage  of 
the  system  is  found  only  in  the  1980  report  of  the  Pipe  Line  Testing 
Company  (5).   Since  this  report  has  only  rough  estimates  of  the  leakage 
flows,  and  since  the  residential  and  commercial  consumption  of  100  gallons 
per  capita  per  day  may  include  leakage  flows,  no  attempt  was  made  to  include 
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Figure  5 

User  Demand  Areas 

o-  Demand  Node 


16 


Nodes 

1-3 
1-4 

3-4 

4-6 
4-6 

6-7 
6-9 
7-8 

8-9 

8-24 

9-11 

11-18 

17-18 

17-19 

18-22 

18-24 

19-20 

19-30 

20-21 

20-30 

21-22 

21-31 

22-23 

23-24 

23-40 
24-25 
25-26 
25-46 
26-47 
26-62 
30-31 
30-33 
31-32 
32-33 
32-40 
33-34 
33-41 
33-34 
34-35 
35-36 
35-38 
35-39 
36-38 
36-39 
39-42 
39-44 
40-46 
40-41 
41-42 
41-50 


Table  3.   Boundary 

Addresses  0 

Boundary 

Addresses 

310-306 

Lexington 

206-205 

Russell 

234-230 

Lexington 

1-0 

Russell 

190-193 

Lexington 

65-59 

Cambridge 

145-144 

Lexington 

162-160 

Bedford 

72-71 

Burlington 

36-34 

Bedford 

1-0 

Lexington 

94-92 

Pleasant 

82-80 

Arlington 

94-96 

Main 

44-42 

Pleasant 

30-31 

N.  Warren 

201-202 

Main 

62-64 

Fowle 

280-281 

Main 

17-18 

Green 

375-380 

Main 

50-55 

Montvale 

438-439 

Main 

529-535 

Main 

75-76 

Salem 

616-620 

Main 

672-673 

Main 

16-17 

Clinton 

41-43 

Eaton 

732-733 

Main 

27-25 

Garfield 

98-100 

Green 

86-88 

Montvale 

159-160 

Montvale 

32-37 

Bow 

210-21 1 

Montvale 

43-42 

Wood 

144-146 

Green 

263-269 

Montvale 

285-286 

Montvale 

76-83 

Central 

22-20 

Central 

45-21 

Washington 

124-127 

Washington 

67-48 

Pine 

133-134 

Washington 

54-52 

Beach 

173-176 

Salem 

235-236 

Salem 

275-300 

Wildwood 

17 


42-44 
44-51 
46-47 
47-48 
48-49 
48-52 
49-50 
49-54 
50-51 
51-58 

52-53 
52-62 

53-54 
53-55 
53-66 
54-56 
55-56 
55-61 

56-57 
57-58 
57-59 
58-60 
59-60 
61-66 

62-63 
63-64 
63-68 
64-65 
64-66 
65-67 
67-70 
67-68 
67-69 
69-70 


Tablr  3 

(continued ) 

380-388 

Salem 

372-376 

Washington 

52-54 

Mishawum 

79-80 

Mishawum 

130-131 

Mishawum 

35-37 

Forest  Park 

83-82 

Olympia 

192-194 

Mishawum 

70-60 

Olympia 

425-444 

Washington 

95-96 

Forest  Park 

24-22 

Alfred 

203-204 

School 

130-1 3^ 

New  Boston 

80-76 

School 

10-20 

Ryan 

20-19 

New  Industrial 

181-183 

New  Boston 

275-280 

Mishawum 

300-310 

Mishawum 

36-39 

Commerce  Way 

9-20 

Dragon  Court 

Commonwealth-Dragon 

Court 

76-75 

Merrimac 

29-30 

Elm 

59-60 

Elm 
West-Pearl 

900-902 

Main 

1  3-1  ^ 

School 
Nichols-Poole 

972-974 

Main 

19-20 

Webster 

39-41 

Poole 

47-45 

Mountain 
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Table  4.   Node  Data 


Ground 

User  Demand-1984 

Location 

Elevation 

(Ft) 

(GPM) 

Node  No 

Samoset  and  Seneca 

161 

181 

1 

Zion  Hill  Reservoir 

302 

0 

2 

Waltham  and  Maura 

235 

90 

3 

Russell  and  Cambridge 

(High  Side) 

64 

106 

4 

Whispering  Hill 

Reservoir 

282 

0 

5 

Russell  and  Cambridge 

(Low  side) 

62 

134 

6 

Bedford  and  Cambridge 

120 

109 

7 

Bedford  and  S.Bedford 

73 

118 

8 

Lexington  and  Garden 

53 

101 

9 

Well  E 

43 

0 

10 

Pleasant  and  Ellis 

81 

73 

11 

Well  D 

56 

12 

Near  Well  F 

47 

0 

13 

Well  F 

46 

0 

14 

Well  C  junction 

60 

0 

15 

Horn  Pond  Reservoir 

239 

0 

16 

Wells  k    ,   B,  C2  47 

Pleasant  and  Arlington  99 

Main  and  Fowle  88 

Main  and  Green  96 

Main  and  Montvale  102 

Main  and  Pleasant  101 

Main  and  Salem  106 

Main  and  Kilby  112 

Main  and  Clinton  120 

Main  and  Wyman  1 10 

North  Warren  and  Ellis  84 

Bedford  and  Winn  126 

Rag  Rock  Reservoir  210 

Green  and  Garfield  124 

Montvale  and  Auburn  127 

Montvale  and  Bow  145 

Montvale  and  Wood  69 

Montvale  and  Green  75 

Montvale  and  Central  62 
Montvale  and  Washington    60 

MDC  Connection  277 
Washington  and  D  Street    56 

Washington  and  Pine  70 

Salem  and  Bow  1  39 

Salem  and  Wood  90 

Salem  and  Pine  77 

Salem  and  Railroad  47 

Cedar  and  Washington  87 

Wells  G  and  H  44 


172 
134 
55 
37 
27 
94 

183 
19 

142 
0 
0 
0 

193 
114 

71 
184 
48 
29 
52 

0 

5 

163 

96 

55 

32 

0 
98 

0 


17 

18 
19 
20 
21 
22 

23 
24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 
44 
45 
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Mishawum  and  Beach  137  40  46 

Mishawum  and  Eaton  144  18  47 

Mishawum  and  Forest  Park  130  51  48 

Mishawum  and  Olympia  81  20  49 

Olympia  and  Wildwood  56  28  50 

Olympia  and  Washington  95  58  51 

Forest  Park  and  Alfred  101  28  52 

Forest  Park  and  School  98  24  53 

Mishawum  and  School  91  13  54 
New  Boston  and  New 

Industrial  96  59  55 

New  Industrial  and  Ryan  81  46  56 


Table  k 

(continued ) 

137 

40 

144 

18 

130 

51 

81 

20 

56 

28 

95 

58 

101 

28 

98 

24 

91 

13 

96 

59 

81 

46 

62 

46 

105 

30 

60 

27 

120 

28 

67 

17 

101 

51 

116 

28 

112 

21 

111 

28 

108 

57 

116 

52 

118 

155 

109 

34 

100 

77 

Mishawum  and  Commerce 

Way  62  46  57 

Mishawum  and  Washington  105  30  58 

Commerce  and  Commonwealth  60  27  59 

Dragon  and  Derry  120  28  60 

New  Boston  and  Merrimac  67  17  61 

Main  and  Alfred  101  51  62 

Elm  and  West  116  28  63 

Elm  and  School  112  21  64 

Main  and  Nichols  111  28  65 

School  and  Merrimac  108  57  66 

Poole  and  Webster  116  52  67 

Webster  and  Pearl  118  155  68 

Winter  and  Mountain  109  34  69 

Main  and  Mountain  100  77  70 


20 


the  spatial  distribution  of  leakage  flows  in  the  user  demand  distribution. 
Instead,  a  "demand  adjustment"  factor  will  be  used  in  the  computer  model 
calibration. 

In  order  to  calculate  the  total  head  at  each  node  the  elevation  of  the 
pipe  is  needed.   The  U.  S.  Geological  Survey  maps  of  Woburn  provide  ten  foot 
elevation  contours.   A  detailed  set  of  topographic  maps  by  Lockwood, 
Kessler,  and  Bartlett  provide  two  foot  elevation  contours.   The  ground 
elevation  for  each  node  was  obtained  from  these  maps  and  is  presented  in 
Table  H.      Figure  6  presents  a  rough  topologic  map  of  Woburn. 

There  are  four  reservoirs  in  the  Woburn  water  distribution  system. 
Their  locations  are  shown  in  Figure  1 .   The  principal  hydraulic  feature  of 
the  reservoirs  is  the  elevation  of  their  water  surface.   Table  5  shows  the 
minimum  and  maximum  elevations,  the  usual  operating  range  of  elevations,  and 
the  storage  volume  of  the  reservoirs. 

While  the  water  pipes  of  Woburn  are  full  of  water  except  for  the  rare 
occasions  when  there  is  a  local  system  failure,  the  storage  reservoirs  have 
regular  emptying  and  filling  cycles.   Typical  cycles  for  the  three 
reservoirs  with  continuously  operating  chart  recorders  are  shown  in  Figures 
7,8,  and  9.   The  Whispering  Hill  reservoir  does  not  have  a  water  level 
measuring  system  in  operation. 

The  cyclic  variation  of  the  reservoir  water  elevation  is  due  to  the 
nature  of  the  user  demand.   The  people  of  Woburn  use  less  water  at  night 
since  most  of  them  are  asleep.   They  use  more  water  in  the  morning  when  they 
bathe  and  wash  clothes.   Industrial  use  also  is  heaviest  in  the  daytime 
since  few  factories  work  three  shifts,  around  the  clock.   People  also  follow 
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50  ft.  Contours 
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Figure  6 

Topographic  Map 
Elevation  in  Feet 
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Table   5.      Reservoir    Data 


Operating        Maximum 
Range  Volume 

Reservoir      Node   No.        Max.    Elev.(Ft)      Min.    Elev.(Ft)        High      Low  (Gal.) 


Zion   Hill 


315 


301  .5 


313        310.5      0.25x10 


6 


Whispering 
Hill 


MDC 


5 


Rag   Rook  29 


Horn   Pond  16 


37 


315 
255 
266 


282.5 

209 

239 


No  data 


4.0x10 


251        2HH  4.4x10 


263        258  6.5x10 


277        277 


The  MDC   connection    is    a  pressure   reducing   valve   set   at   94    psi    or   217   ft    at   a 
ground   elevation  of  60    ft.      The  setting   can  be   varied   by   the  operator. 
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Figure  7 

Zion  Hill  Reservoir  Chart 
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Figure  8 

Rag  Rock  Reservoir  Chart 
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Figure  9 

Horn  Pond  Reservoir  Chart 
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weekly  and  annual  cycles.   The  work  week  is  usually  Monday  to  Friday.   Crops 
and  gardens  and  swimming  pools  demand  extra  water  during  the  summer. 

As  can  be  seen  from  Figure  8,  the  weekly  cycle  of  filling  and  emptying 
the  Rag  Rock  water  tank  is  weak  while  the  daily  cycle  is  strong.   In  the 
night,  say  midnight  to  6  AM,  the  water  level  rises  three  or  four  feet  with  a 
corresponding  storage  volume  increase  of  300,000  or  ^00,000  gallons.   The 
filling  flow  rate  is  thus  1.2  to  1  .6  million  gallons  per  day.   The  daytime 
emptying  flow  rate  is  similar  in  the  morning  and  less  in  the  evening.   In 
the  summertime  the  daily  cycle  is  stronger. 

The  daily  filling  and  emptying  cycle  of  the  Rag  Rock  tank  is  the 
principal  indicator  of  the  peaking  factor  of  user  demand.   An  examination  of 
many  weekly  Rag  Rock  charts  shows  that  the  maximum  and  minimum  demands  are 
roughly  20  percent  above  and  below  the  daily  average  demand.   Therefore  the 
corresponding  peaking  factors  for  the  daily  cycle  are  1.2  and  0.8. 

The  other  storage  reservoirs  of  Woburn  are  not  as  affected  by  these 
cycles.   The  Zion  Hill  and  Whispering  Hill  reservoirs  are  maintained  at  a 
nearly  constant  level  by  the  Shaker  Glen  booster  pumps  shown  in  Figure  6. 
Since  the  reservoirs  are  located  in  the  highest  part  of  Woburn,  these  pumps 
are  needed  to  increase  the  total  head  to  the  level  of  these  reservoirs  for 
the  Shaker  Glen  consumers  they  serve.   When  the  reservoir  levels  drop  only 
slightly  (1.5  ft)  from  the  normal  full  levels,  the  booster  pump  turns  on. 
Thus  the  Zion  Hill  recorder's  chart,  Figure  7,  shows  only  this  small  level 
change. 

The  Horn  Pond  reservoir  is  so  close  to  Wells  A?,    B,  C?,  and  D  that  its 

elevation  is  closely  set  by  the  pumping  head  and  throttle  valves  of  these 
wells.   Figure  9  shows  the  pressure  of  the  water  in  the  pipe  near  Well  Cp 
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and  thus  indicates  the  total  head  in  che  Horn  Pond  reservoir.   Although 
there  is  a  small  variation  when  Well  C      is  turned  on  and  off,  this  pressure 

is  a  good  indicator  of  the  water  level  in  the  Horn  pond  reservoir.   As 
Figure  9  shows,  it  is  nearly  constant. 

Besides  the  four  actual  reservoirs  in  Woburn  there  is  also  another 
hydraulic  element  that  will  be  treated  as  if  it  were  a  reservoir.   The  MDC 
water  supply  connection  to  Woburn  occurs  at  Montvale  Avenue  and  Washington 
Street,  where  a  pressure  regulating  valve  changes  the  high  pressure  of  the 
^8  inch  pipe  of  the  MDC  to  the  lower,  regulated  pressure  of  Woburn's  12  inch 
pipe.   The  valve  holds  the  pressure  steady  and  supplies  whatever  water  flow 
is  needed  to  do  so,  thus  acting  as  if  it  were  a  reservoir  supplying  water  at 
a  constant  total  head. 

There  are  eight  well  pumps  in  Woburn  and  two  principal  booster  pumps. 
The  well  pumps  supply  the  energy  (total  head)  required  to  raise  the  ground 
water  from  its  level  in  the  wells  to  the  reservoirs  and  users  of  the  city. 
The  pumps  are  described  by  their  characteristic  curves  of  pump  head  vs.  flow 
rate  and  their  operating  points  on  those  curves.   The  characteristic  curves 
for  the  Layne  &  Bowler  pumps  (Z)  of  wells  E,  F,  G,  and  H  are  shown  in 
Figures  10,  11,  12,  and  13  respectively.   Although  there  is  some  variation 
they  typically  have  an  operating  point  at  a  head  of  265  feet.   The  other 
four  pumps,  for  wells  A  ,  B,  C?   and  D  cannot  be  described  by  their 

characteristic  curves  since  they  are  controlled  by  throttle  valves.   These 
valves  are  used  by  the  Woburn  pump  operators  to  control  those  pumps  and  get 
the  flow  rate  that  the  operators  desire.   Therefore  they  are  simply 
characterized  by  their  flowrate.   Under  normal  operating  conditions  the 
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Figure  10a 
Well   E  1967 
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Figure  11a 
Well  F  1960 
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Figure  lib 
Well  F  1973 
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Figure   12 
Well  G  1978 
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Figure  13 
Well  H  1977 
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positions  of  the  throttle  valves  are  unchanged;  wells  A„  ,  B,  and  D  are  run 

continuously;  and  well  C?  is  run  intermittently  to  provide  the  total  flow 

needed  to  maintain  the  Horn  Pond  reservoir  level  near  262  feet.   The  two 
booster  pumps  had  well  defined  design  operating  conditions  but  poorly 
defined  characteristic  curves.   The  only  curve  available  was  on  a  pump 
selection  chart,  Figure  14. 

The  pumps  need  more  maintenance  in  order  to  maintain  their  operating 
characteristics  than  do  the  other  elements  of  the  water  system.   Well  E  was 
repaired  in  1967  and  1977,  with  a  significant  improvement  from  350  gpm  to 
670  gpm  being  documented  for  the  1977  repair.   Well  F  was  considered  a  350 
gpm  -  240  ft  head  pump  in  1960  but  was  repaired  and  modified  in  1973  to 
serve  as  a  300  gpm  -  328  ft  head  pump.   Wells  G  and  H  were  repaired  in  1978 
and  1977  respectively,  with  an  improvement  in  flow  for  well  H  from  150  to 
400  gpm  being  recorded.   Thus  the  well  characteristics  shown  in  Figures  10 
to  14  should  be  considered  the  "new"  curves  and  the  effects  of  aging  should 
be  taken  into  account.   This  could  not  be  done  due  to  the  lack  of  data  on 
the  history  of  the  pump  curves.   Table  6  summarizes  the  basic 
characteristics  and  repair  histories  of  the  four  well  pumps  without  throttle 
valves. 

The  pumps  of  Woburn  are  all  metered  and  so  these  meters  provide  the 
data  needed  to  document  the  volume  of  water  pumped  each  minute  and  the 
operating  procedure  of  the  water  system,  as  well  as  the  aging  of  the  pumps. 
Most  of  the  pumps  operate  continuously  and  steadily,  the  same  flow  rate  hour 
after  hour,  day  after  day.   The  Shaker  Glen  booster  pumps  do  not,  turning  on 
and  off  regularly  as  determined  by  the  water  level  in  the  Zion  Hill  tank. 


35 


150 

*:  ioo 

_     80 


O 


60 


-J I l_ 

400        800 


Flow  Rate  Q  gpm 


Figure   14 
Booster  pump 
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Table  6.   Pump  Characteristics  and  History 


Pump 

Well  E 


Design  Flow  Rate  (GPM)    Design  Head  (ft) 


Well  F 


Well  G 


Well  H 


830 
600 
600 

350 
300 

800 
800 

400 

400 


263 
264 

240 
328 

265 
265 

265 
264 


Date  Installed 
or  Repaired 

1937  (installed) 
1966 

1977 

1960  (installed) 
1973 

1964  (installed) 
1978 

1964  (installed) 
1977 
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Also,  as  previously  mentioned,  well  C?  does  not  operate  continuously.   It  is 

turned  on  and  off  manually  by  the  pump  station  operator,  whose  decisions  are 
based  on  the  water  levels  in  the  Horn  Pond  reservoir.   Thus  the  booster 
pumps  are  used  to  automatically  maintain  a  close  balance  between  water 
supply  and  water  demand  in  the  Shaker  Glen  high  pressure  zone  while  well  C 

is  controlled  manually  to  provide  the  same  balance  in  the  overall  system. 

The  well-documented  on/off  cycles  of  the  control  pumps  cause  a  wider 
variety  of  flow  patterns  than  those  considered  by  Waldorf  and  Cleary  (6). 
The  effects  of  these  cycles  was  investigated  during  the  model  calibration 
phase  of  this  project  and  will  be  discussed  in  that  section  of  this  report. 

The  annual  cycle  of  water  pumpage  is  shown  in  Table  7  and  Figure  15. 
Since  these  pump  records  involve  much  more  water  than  can  be  stored  in 
Woburn's  reservoirs,  these  water  supply  figures  also  indicate  the  annual 
cycle  of  user  demand.   The  winter  requires  ten  percent  less  water  than  the 
annual  average  when  the  summer  needs  15  percent  more.   This  extra  demand 
during  the  summer  taxed  the  limits  of  the  supply  system's  capacity  and  water 
use  restrictions  were  needed  at  times  of  generally  low  water  levels  in  the 
reservoirs. 

Tables  8  and  9  are  taken  directly  from  the  report  of  Waldorf  and  Cleary 
(6).   They  show  the  percent  of  the  total  well  flow  that  came  from  wells  G 
and  H  and  the  number  of  days  when  these  two  wells  were  pumping.   Realizing 
that  the  G  and  H  flow  percentage  varied  substantially,  Waldorf  and  Cleary 
selected  four  typical  months  for  each  of  the  historical  periods,  1964-1969 
and  1970-1979.   They  then  calculated  the  four  zones  where  water  from  wells  G 
and  H  was  used  during  those  four  typical  months.   This  report  attempts  to 
follow  this  procedure  of  Waldorf  and  Cleary.   Their  pumping  data  for  the 
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Table  7. 

Monthly 

Averag 

;e  Pump 

Rate  in  MGD 

Year 

Jan 

Feb 

Mar 

April 

May 

June 

July 

Aug 

Sept. 

Oct. 

Nov. 

Dec. 

1964 

3.58 

3.93 

4.33 

4.45 

5.67 

6.15 

4.92 

4.87 

4.32 

3.93 

3.59 

3.62 

1965 

3.54 

3.74 

3.80 

3.83 

4.69 

5.30 

5.86 

5.15 

3.89 

3.79 

3-67 

3.42 

1966 

3.51 

3.55 

3.65 

4.28 

4.18 

4.88 

.6.08 

5.58 

4.45 

4.12 

4.10 

3.75 

1967 

3.85 

4.27 

3.82 

4.21 

4.28 

5.43 

4.55 

5.02 

5.14 

4.94 

4.84 

4.80 

1968 

4.80 

4.65 

4.71 

5.16 

5.32 

5.25 

6.13 

6.26 

5.68 

5.25 

4.88 

4.76 

1969 

4.94 

5.00 

5.00 

5.00 

5.27 

6.49 

6.30 

6.05 

5.32 

4.71 

4.59 

4.57 

1970 

4.59 

4.60 

4.71 

4.82 

5.16 

6.35 

6.09 

5.96 

4.91 

4.66 

4.38 

4.27 

1971 

4.64 

4.64 

4.57 

4.44 

4.72 

5.92 

6.38 

5.53 

5.08 

4.54 

4,21 

4.20 

1972 

4.51 

4.36 

4.46 

4.56 

5.02 

5.36 

5.55 

5.38 

4.92 

4.53 

4.36 

4.31 

1973 

1974 

4.74 

4.84 

4.84 

4.83 

5.20 

5.58 

5.69 

6.33 

5.20 

5.10 

4.78 

4.36 

1975 

4.83 

4.93 

4.83 

4.76 

5.47 

5.89 

6.35 

5.77 

5.12 

4.93 

5.08 

5.07 

1976 

4.92 

5.36 

5.20 

5.53 

5.51 

7.00 

6.07 

5.68 

5.51 

5.51 

5.09 

4.93 

1977 

5.93 

5.70 

5.63 

5.29 

5.69 

5.61 

6.00 

5.68 

5.61 

5.27 

5.21 

5.30 

1978 

5.40 

5.43 

5.48 

5.68 

5.99 

6.81 

6.98 

6.18 

6.04 

5.588 

5.60 

5.94 

1979 

'5.94 

6.24 

'6.46 

6.05 

6.02 

6.98 

Avg. 

4.65 

4.75 

4.77 

4.86 

5.21 

5.93 

5.93 

5.76 

5.09 

4.78 

4.60 

4.52 
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Figure  15 

Annual  Cycle 
1964  - 1979 
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Ta  b  le 

8.   Wells  G  an 
X 

d  H  Pumping 
G+H 

His  tory  : 

°lHxl00 
X 

1964-1969 

(6) 

Average 

Rate,  Wells 

Daily  Supply 

Pumped 

%  G+H 

Zone 

//Days  Wells 

Year 

Month 

in  MGD 

in  MGD 

by  Volume 

Category 

G+H  Pumped 

(all  wells) 

1964 

Oct 

3.93 

0.715 

18.2 

2 

31 

Nov 

3.59 

0.705 

19.6 

2 

24 

Dec 

3.62 

0.689 

19.0 

2 

28 

1965 

Dec 

3.42 

1.202 

35.1 

4 

30 

1966 

Jan 

3.51 

1.250 

35.6 

4 

31 

Feb 

3.55 

1.150 

32.4 

4 

28 

Mar 

3.65 

1.050 

28.8 

3 

31 

Apr 

4.28 

0.999 

23.3 

2 

30 

May 

4.18 

0.975 

23.3 

2 

30 

Jun 

4.88 

1.000 

20.5 

2 

30 

Jul 

6.08 

1.000 

16.4 

1 

31 

Aug 

5.58 

1.000 

17.9 

2 

31 

Sep 

4.45 

1.000 

22.5 

2 

30 

Oct 

4.12 

1.000 

24.3 

2 

31 

Nov 

4.10 

0.852 

20.8 

2 

20 

Dec 

3.75 

0.822 

21.9 

2 

31 

1967 

Jan 

3.85 

0.822 

21.4 

2 

31 

Feb 

4.27 

1.000 

23.4 

2 

28 

Mar 

3.82 

1.000 

26.2 

3 

12 

Jul 

4.55 

0.576 

12.7 

1 

3 

Aug 

5.02 

0.576 

11.5 

1 

8  H  Alone 

1.728 

34.4 

4 

20  H+G  Together 

Sep 

5.14 

1.152 

22.4 

2 

1  G  Alone 

1.728 

33.6 

4 

26  G+H  Together 

Oct 

4.94 

1.728 

35.0 

4 

31 

Nov 

4.84 

1.438 

29.7 

4 

30 

Dec 

4..  80 

1.491 

31.1 

4 

5 

1968 

Apr 

5.16 

1.152 

22.3 

2 

19 

May 

5.32 

1.152 

21.7 

2 

31 

Jun 

5.25 

1.152 

21.9 

2 

30 

Jul 

6.13 

1.152 

18.8 

2 

31 

Aug 

6.26 

1.152 

18.4 

2 

31 

Sep 

5.68 

1.148 

20.2 

2 

30 

Oct 

5.25 

1.152 

21.9 

2 

31 

Nov 

4.88 

0.954 

19.5 

2 

30 

Dec 

4.76 

0.948 

19.1 

2 

31 

1969 

Jan 

4.94 

0.938 

19.0 

2 

31 

Feb 

5.00 

1.152 

23.0 

2 

28 

Mar 

5.00 

1.152 

23.0 

2 

31 

Apr 

5.00 

1.152 

23.0 

2 

30 

May 

5.27 

1.152 

21.9 

2 

31 

Jun 

6.49 

1.152 

17.8 

2 

30 

Jul 

6.30 

1.152 

18.3 

2 

31 

Aug 

6.05 

1.152 

19.0 

2 

31 

Sep 

5.32 

1.152 

21.7 

2 

30 
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Summary  of  1964  to  1969 
Pumping  Data 


Total  days  wells  G+H  pumped   -  1,199 
Maximum  G+H  Percent  by  Volume  -  35.6 
Minimum  G+H  Percent  by  Volume  -  11.5 

Zone  Boundaries 

%  G+H  by 
Zone  #       Volume    Representative  Month 


1 

11.5-17.5 

June    '69 

2 

17.5-23.5 

Feb.    f67 

3 

23.5-29.5 

Nov.    '67 

4 

29.5-35.6 

Jan.    '66 
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Table 

9.   Wells  G  and  H  Pumping 

History: 

1970-1979 

(6) 

X 

G+H 

-Y~xl00 

Average 

Rate,  Wells 

Daily  Supply 

Pumped 

%   G+H 

Zone 

//Days  1 

rfells 

Year 

Month 

in  MGD 

in  MGD 

by  Volume 

Category 

G+H  Pump 

ed 

(all  wells) 

1970 

Apr 

A.  82 

0.576 

12.0 

1 

30 

May 

5.16 

0.808 

15.7 

2 

31 

Jun 

6.35 

1.152 

18.1 

2 

30 

Jul 

6.09 

1.135 

18.6 

2 

31 

Aug 

5.96 

1.152 

19.3 

2 

31 

Sep 

4.91 

1.152 

23.5 

3 

30 

Oct 

4.86 

0.576 

11.9 

1 

31 

Nov 

4.88 

0.381 

7.8 

1 

30 

Dec 

4.27 

0.384 

9.0 

1 

31 

1971 

May 

4.72 

0.400 

8.5 

1 

9 

Jul 

6.39 

0.908 

14.2 

1 

29 

Aug 

5.53 

1.147 

20.7 

2 

31 

Sep 

5.09 

1.110 

21.8 

3 

29 

Oct 

4.54 

1.132 

24.9 

3 

26 

Nov 

4.21 

0.907 

21.5 

3 

21 

Dec 

4.20 

0.937 

22.3 

3 

31 

1972 

Jan 

4.51 

0.957 

21.2 

3 

31 

Feb 

4.36 

1.152 

26.4 

3 

29 

Mar 

4.45 

1.126 

25.3 

3 

29 

Apr 

4.56 

1.012 

22.2 

3 

17 

1974 

Aug 

6.33 

1.595 

25.2 

3 

4 

H+G  Together 

1.069 

16.9 

2 

27 

G 

Alone 

Sep 

5.20 

1.545 

29.7 

4 

30 

Oct 

5.10 

1.002 

19.6 

2 

1 

G 

Alone 

1.581 

31.0 

4 

30 

H+G  Together 

Nov 

4.78 

1.008 

21.1 

3 

10 

G 

Alone 

1.553 

32.5 

4 

20 

H 

and  G 

Dec 

4.36 

0.999 

22.9 

3 

8 

G 

Alone 

1.471 

33.7 

4 

5 

H 

and  G 

1975 

May 

5.47 

1.047 

19.1 

2 

15 

Jun 

5.89 

1.368 

23.2 

3 

4 

H 

and  G 

1.037 

17.6 

2 

26 

G 

Alone 

Jul 

6.35 

1.560 

24.6 

3 

31 

Aug 

5.77 

1.008 

17.5 

2 

13 

G 

Alone 

1.582 

27.4 

4 

18 

H 

and  G 

Sep 

5.12 

1.533 

29.9 

4 

10 

H 

and  G 

0.996 

19.5 

2 

20 

G 

Alone 

Oct 

4.93 

0.977 

19.8 

2 

20 

Nov 

5.08 

1.004 

19.8 

2 

13 

Dec 

5.07 

0.984 

19.4 

2 

31 

1976 

Jan 

4.92 

1.052 

21.4 

3 

31 

Feb 

5.36 

1.022 

19.1 

2 

29 

Mar 

5.20 

0.971 

18.7 

2 

18 

Jun 

7.00 

1.531 

21.9 

3 

22 

Jul 

6.87 

1.608 

23.4 

3 

31 

Aug 

6.07 

0.576 

9.5 

1 

8 

H 

Alone 

1.512 

24.9 

3 

23 

H 

and  G 

Sep 

5.68 

0.576 

10.1 

1 

30 

Oct 

5.51 

0.576 

10.5 

1 

30 
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Table  9  (continued) 


Year 


1977 


1978 


1979 


Average 
Daily  Supply 
Month    in  MGD 

(all  wells) 
Jan  5.93 
Feb  5.70 
Mar  5.63 
Apr  5.29 
May      5.69 

Jun  5 . 81 


Jul 

6.00 

Aug 

5.68 

Sep 

5.61 

Oct 

5.27 

Nov 

5.21 

Dec 

5.30 

Jan 

5.40 

Feb 

5.43 

Mar 

5.48 

Apr 

5.68 

May 

5.99 

Jun 

6.81 

Jul 

6.98 

Aug 

6.18 

Sep 

6.04 

Oct 

5.58 

Nov 

5.60 

Dec 

5.94 

Jan 

5.94 

Feb 

6.31 

Mar 

6.46 

Apr 

6.05 

May 

6.42 

G+H 
Rate,   Wells 
Pumped 
in  MGD 

1.044 

1.029 

1.010 

1.007 

1.639 

1.029 

1.496 

1.002 

0.991 

0.976 

0.963 

0.953 

0.947 

0.970 

0.980 

0.978 

0.927 

1.423 

1.662 

1.627 

1.629 

1.635 

1.612 

1.618 

1.617 

1.624 

1.57 

1.723 

1.575 

1.548 

1.600 


£±Kxl00 
X 

%  G+H 
by  Volume 

17.6 

18.1 

17.9 

19.0 

28.8 

18.1 

25.7 

17.2 

16.5 

17.2 

17.2 

18.1 

18.2 

18.3 

18.1 

18.0 

16.9 

25.1 

27.7 

23.9 

23.3 

26.5 

26.7 

29.0 

28.9 

27.3 

26.4 

27.3 

24.4 

25.6 

24.9 


Zone     //Days  Wells 
Category  G+H  Pumped 


2 

2 

2 

2 

4 

2 

3 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

4 

3 

3 

3 

3 

4- 

4 

4 

3 

3 

3 

3 

3 


10 
28 

29 
30 

12  H+G 


19 
3 
27 
31 
31 
30 
31 
30 
31 
31 
28 
31 
30 
31 
30 
31 
31 
30 
31 
30 
31 
31 
28 
31 
30 
21 


G  only 
H+G 
G  only 


Summary  of  1970  to  1979 
Pumping  Data 


Total  days  wells  G+H  pumped   -  1,963 
Maximun  G+H  percent  by  volume  -  33.7 
Minimum  G+H  percent  by  volume  -  7.8 

Zone  Boundaries 


Zone  # 

%  G+H  by 
Volume 

Representative  Month 

1 
2 
3 
4 

7.8-14.3 
14.3-20.8 
20.8-27.3 
27.3-33.7 

July  '71 
Aug.  '71 
Feb.  '79 
Dec.  '74 
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eight  typical  months  are  given  in  Table  10.   However,  the  pumping  history 
given  by  Tables  8  and  9  was  found  to  be  erroneous  and  so  a  significant 
change  in  the  analysis  of  the  pumping  history  of  wells  G  and  H  was  required. 
This  new  analysis  is  presented  in  Section  V. 

In  addition  to  the  data  available  from  the  City  of  Woburn,  substantial 
data  was  also  supplied  by  the  ISO  Commercial  Risk  Services  Company.   This 
company  gathers  data  on  water  distribution  systems  in  order  to  determine  a 
fire  insurance  classification  which  is  used  to  develop  fire  insurance  rates. 

The  ISO  measurement  procedure  uses  at  least  two  nearby  fire  hydrants. 
Usually  a  pressure  gauge  mounted  on  a  pipe  plug  is  attached  to  the  hydrant 
nearer  to  the  source  of  water,  the  hydrant  valve  is  opened  and  the  water 
pressure  is  recorded.   Then  the  valve  of  the  farther  hydrant  is  opened  and  a 
jet  of  water  issues  forth  from  this  hydrant.   The  velocity  of  the  jet  is 
measured  using  a  special  Pitot  tube  and  pressure  gauge.   At  the  same  time, 
with  the  jet  flowing  from  the  farther  hydrant,  the  pressure  gauge  on  the 
nearer  hydrant  is  read  again. 

An  extensive  set  of  ISO  measurements  was  recorded  on  November  2  and  3» 
1983.   This  data  is  presented  in  Table  11  and  is  used  to  calibrate  the 
computer  model.   Not  all  of  this  1983  ISO  data  is  useful.   The  data  from 
"dead  end"  pipes  and  other  pipes  eliminated  from  the  network  in  the 
simplification  process  cannot  be  directly  compared  with  the  model.   However 
the  total  head  of  all  the  ISO  tests  can  be  calculated  and  used  in  a  general 
comparison  with  the  model.   The  elevation  of  the  hydrants  is  the  only 
additional  data  needed  to  find  these  total  heads  since  the  velocity  head  is 
negligible. 
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Table  10.   Typical  Months  of  Pumping  in  MGD  (6) 

1964-1969 

Zone  I        II       III      IV 

Month 
Total  Flow 
%   G  &  H 

Wells  A2,  B,  C2 

Well  D 
Well  E 
Well  F 
Wells  G  &  H 


6/69 

2/67 

1 1/67 

1/66 

6.72 

4.25 

4.86 

3.55 

17.1 

23.5 

29.6 

35.2 

3.70 

2.92 

1  .95 

2.30 

0.90 

0. 

0.94 

0. 

0.65 

0. 

0.53 

0. 

0.32 

0.33 

0. 

0. 

1.15 

1  .00 

1  .43 

1.25 

1970- 

-1979 

Zone  I        II       III      IV 

Month 
Total  Flow 
%G   &  H 

Wells  Ap,  B,  C2 

Well  D 
Well  E 
Well  F 
Wells  G  &  H 


7/71 

8/71 

2/79 

12/74 

6.41 

5.49 

6.31 

4.41 

14.2 

20.9 

27.3 

33.4 

3.56 

3.40 

3-37 

1  .91 

0.99 

0. 

0.73 

0.43 

0.67 

0.68 

0.04 

0.60 

0.28 

0.26 

0.45 

0. 

0.91 

1  .15 

1  .72 

1  .47 
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Table  11.   1983  ISO  Data 


Location 


Hudson  and  Arlington 
Cranes  Court  off  Main 
Main  and  Montvale 
Jefferson  and  Prospect 
Garfield  and  Sonar 
Micro  and  Holton 
Henshaw  and  Erie 
Montvale  and  Campbell 
Unicorn  Park  and  Montvale 
Montvale  Road  and  Dale 
Salem  St.  and  Salem  Ave. 
Cedar  and  Washington 
Dragon  and  Derry 
Commerce  and  Mishawum 
Commerce  and  Atlantic 
School  and  New  Boston 
New  Boston  and  Sixth 
New  Boston  @  end 
Foster  and  Jones 
Main  and  Mountain 
Webster  St.  and  Webster  Ave, 
Westgate  Drive  and  Pearl 
Middlesex  Canal  and  Main 
Main  and  Clinton 
Continental  Court  and  Salem 
Middle  and  Dorothy 
Mishawum  and  Olympia 
Bedford  and  Houghton 
Bedford  and  Cambridge 
Locust  and  Willow 
Totman  and  Lexington 
Frances  and  Main 
Cambridge  and  Country  Club 
Lexington  and  Cambridge 

(High  Si^e) 
Carol  and  Anna 
Russell  and  Silvermine 


Hydrant 
Flow  (GPM) 

2550 
2125 
1965 
1670 
630 
1800 
1575 
1470 
2050 

1625 
1130 
2025 
1395 
1075 
1545 
1545 
1610 
1290 

650 
1010 
1625 

780 
2075 
2550 
1965 
1475 
1545 
2225 
1670 
1230 

380 

745 

970 

2550 

1280 

925 


Pressure  (psi) 

Flow 

No  Flow 

76 

83 

62 

74 

59 

62 

40 

64 

21 

55 

59 

82 

60 

94 

51 

54 

40 

94 

64 

75 

10 

70 

58 

71 

34 

54 

73 

81 

58 

76 

55 

69 

64 

74 

36 

74 

28 

64 

25 

58 

26 

52 

21 

55 

46 

57 

50 

54 

36 

44 

30 

51 

67 

69 

42 

60 

18 

55 

44 

78 

12 

77 

42 

52 

47 

64 

100 

105 

38 

60 

12 

49 

Node  No, 


21 


44 
60 
57 
59 
53 

61 

70 
67 


25 

40 

49 
7 

24 


4 

1 
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While  this  ISO  pressure  data  is  very  good  for  model  calibration,  there 
is  one  problem  that  was  "discovered"  during  the  validation  process.   The 
hydrants  are  not  neatly  located  in  pairs  near  the  principal  pipe  junctions 
of  the  computer  model.   Thus  an  ISO  location  like  Main  Street  and  Montvale 
Avenue  may  refer  to  a  variety  of  hydrants  near  that  pipe  junction.   If  the 
ground  is  fairly  flat  near  the  junction  this  causes  no  problem.   When  the 
various  hydrants  are  at  substantially  different  elevations,  then  the  total 
head  associated  with  the  ISO  measurement  may  be  in  error  due  to  the  use  of  a 
pressure  from  one  hydrant  and  an  elevation  from  another.   Since  the  exact 
hydrants  tested  are  unknown,  the  ground  elevation  of  the  pipe  junction  was 
used  in  figuring  total  heads  and  the  possible  error  noted. 

Some  Woburn  system  data  that  complements  the  1983  ISO  data  is  presented 
in  Table  12.   The  reservoir  and  pump  charts  and  other  sources  of  this 
complementary  data  are  given  in  Appendix  III. 

In  addition  to  the  1983  ISO  data  used  for  the  computer  model 
calibration,  there  is  also  available  data  for  the  1964-1979  period  when 
wells  G  and  H  were  in  operation.   This  data  is  also  included  in  Appendix  IV. 
III.   Computer  Model 

The  analysis  of  a  water  distribution  system  is  done  by  applying  the 
basic  principles  of  the  conservation  of  mass  and  energy  to  the  water  flowing 
through  the  system.   The  mass  conservation  principal  is  applied  by  assuming 
that  the  pipes  do  not  leak  along  their  length  and  that  the  users  take  water 
out  of  the  system  only  at  pipe  junctions.   Thus  the  net  flow  rate  into  each 
pipe  junction  must  equal  the  user  demand  of  the  junction.   UsuaLly  the  user 
demand  is  known  and  the  flow  rates,  Q,  in  the  pipes  are  the  unknowns  to  be 
determined  by  the  analysis.   The  energy  conservation  principal  is  expressed 
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Table  12.   Water  System  Data  for  November  2  &  3,  1983 


Reservoirs 

Horn  Pond  Elevation:  259  to  260  ft 

Rag  Rock  Elevation:  249  to  245  ft 

Outflow:  1250  to  0.  gpm 

Zion  Hill  Elevation:  311  to  313  ft 

MDC  Pressure:  94  psi 

Pumps 

Well  A0  Q  =  1085  gpm    (steady  flow  -  2  day  avg.) 

Well  B  Q  =  173  gpm     (monthly  avg.) 

Well  C2  Q  =  950  gpm     (on  from  2:15  to  10:50  PM, 

11/2) 
(on  from  8:00  AM  to  2:00  PM, 

11/3) 
Well  D  Q  =  400  gpm     (steady  flow  -  4  day  avg.) 

Well  E  Q  =  454  gpm     (monthly  avg.) 

Well  F  Q  =  262  gpm     (monthly  avg.) 

Shaker  Glen  Booster  "on" 
MDC  Q  =  1604  gpm    (monthly  avg.) 
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in  terms  of  height,  usually  referred  to  as  "head".   Potential  energy  is 

described  by  the  elevation,  Z,  of  the  centerline  of  the  pipe  above  mean  sea 

level.   The  kinetic  energy  per  unit  weight  is  given  by  the  square  of  the 

2 
flow  velocity  divided  by  twice  the  acceleration  of  gravity,  V  /2g.   The 

pressure  energy  per  unit  weight  is  given  by  P/Y,  where  P  is  the  pressure  and 

Y  is  the  specific  weight  of  the  water.   All  these  three  kinds  of  energy  are 

combined  to  form,  H,  the  "total  head"  or  mechanical  energy  per  unit  weight 

by  the  equation, 

H  =  Z  +  V2/2g  +  P/Y  (1) 

The  energy  principal  is  applied  to  pipes  by  requiring  that  the  energy  lost 
to  friction,  h  ,  is  the  same  as  the  energy  difference  between  the  junctions 

at  the  ends  of  the  pipe.   Since  the  flow  in  pipes  in  series  is  easily 
determined,  only  those  pipes  that  parallel  other  pipes,  forming  loops,  need 
this  energy  analysis.   A  corresponding  energy  analysis  is  applied  to  pumps 
and  reservoirs  and  valves  also. 

The  result  is  a  large  system  of  non-linear  equations  which  are  solved 
to  determine  the  flow  rates  in  each  pipe.   More  detail  on  the  theory  of 
water  distribution  system  analysis  can  be  found  in  elementary  fluid 
mechanics  textbooks  and  in  Analysis  of  Flow  in  Pipe  Networks  by  Roland  W. 
Jeppson  (2). 

The  computer  program  for  the  pipe  network  analysis  of  this  study  is 
called  NETWK.   It  was  developed  by  Professor  Roland  W.  Jeppson  of  Utah  State 
University  in  association  with  the  consulting  engineering  firm,  CH?M  Hill  of 

Corvallis,  Oregon.   It  is  available  for  teaching  and  research  purposes  at 
the  University  of  Massachusetts  Computing  Center  and  was  used  there. 
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As  stored  in  the  University  Computer  Center,  the  program  can  handle  up 
to  500  pipes,  230  nodes,  20  reservoirs,  20  pumps,  20  pressure  reducing 
valves  and  90  minoi — loss  devices.   The  data  requirements  for  the  pipes  are 
the  length,  diameter,  and  friction  coefficient  of  each  pipe,  and  the  nodes 
to  which  the  pipes  are  connected.   The  user  demand  and  elevation  of  each 
node  must  be  given.   Reservoirs  are  described  by  the  elevation  of  their 
water  level.   Three  points  from  the  characteristic  curve  of  each  pump  are 
required,  as  well  as  the  elevation  of  the  water  on  the  suction  side  of  the 
pump.   A  more  detailed  list  of  the  input  data  requirements  of  the  NETWK 
program  can  be  found  in  the  Users'  Manual  (3). 

The  first  version  of  the  input  data  for  the  computer  model  is  given  in 
Table  13.   This  input  file  began  the  model  calibration  process.   The  file 
name  is  W0BURNA,  under  the  account  name  ADWE000.   The  pipes  are  described  by 
the  pipe  number,  the  two  pipe  junction  or  node  numbers,  and  the  pipe  length 
(feet),  diameter  (inches),  and  roughness  coefficient.   The  nodes  or  pipe 
junctions  are  given  by  their  node  number,  user  demand  (gallons  per  minute), 
and  ground  elevation  (feet)  above  mean  sea  level.   The  flows  for  wells  A  , 

B,  and  C  and  for  well  D  are  included  as  negative  user  demands  at  nodes  17 

and  12  respectively.   The  reservoirs  are  listed  by  their  node  numbers  and 
water  surface  elevations  (feet).   The  well  pumps  are  described  by  their  node 
number,  three  sets  of  flow  rate-pump  head  (GPM,  feet)  points  from  their 
characteristic  curves,  and  the  ground  water  level  (feet)  in  the  wells.   The 
booster  pump  is  similar  to  the  wells,  except  that  it  is  identified  by  its 
pipe  number  and  its  water  level  is  zero.   The  check  valve  in  the  booster 
pump  pipe  is  simply  identified  by  its  pipe  number. 
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Table  13.      Computer    Model    Input    Data    (see    P.    51) 

p* 

WATER  DISTRIBUTION  IN  UOBURNn  MA 

MODEL  CALIBRATION?  NOVEMBER  ZvA*     1984 
/* 
PIPES 

1  4   1  4200  6  30  / 

2  3   1  5200  10  SO  / 

3  4   3  3700  12  30  / 

4  3   2  700  10  80  / 

5  4   5  4600  12  30  / 

6  6   4  400  12  30  / 

7  6  7  4800  6  30  / 
3  3  7  2300  10  30  / 
9  9   6  3000  6  30  / 

10  9   3  4500  12  30  / 

11  10   9  800  12  30  / 

12  9  11  2600  6  SO  / 

13  24   8  7000  10  30  / 

14  12   6  7200  12  30  / 

15  12  13  5700  16  80  / 

16  13  11  3700  16  80  / 

17  14  13  100  10  30  / 
13  12  15  300  14  30  / 

19  15  16  500  16  80  / 

20  15  17  200  16  30  / 

21  17  18  6000  16  30  / 

22  11  18  1300  16  80  / 

23  17  19  4700  16  30  / 

24  19  20  900  16  80  / 

25  20  21  2000  16  30  / 
">c  21  22  300  16  30  / 


27  18  22  1300  16  80  / 

28  22  23  900  16  30  / 

29  23  24  1400  16  30  / 

30  24  25  1900  16  80  / 

31  25  26  700  16  80  / 

32  11  27  1700  16  80  / 

33  13  27  1100  10  SO 

34  22  27  2100  16  30 

35  27  28  1750  20  80 

36  23  29  1400  24  30  / 

37  28  26  3800  16  80  / 

38  19  30  1500   3  30  / 

39  20  30  1700  10  30  / 

40  21  31  1100  12  30 

41  30  31  3400   8  SO 


*j  v  *.>  .i.      ».>  ■■*  \y  \/         i_>    o  \;     / 
4 *.  Jl   J .-..       .1.  >  U  O   I  v_     O  O   / 


43  33  32    2200  1.2  80 

44  34  ZZ  1000  12  SO 

45  34  30    4500  10  80 

47  Z6    35     800  12  80 

48  37  Z6  100  12 

52 


o  r. 


Table    13    (continued) 
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35 

38 

3300 

80 

50 

36 

38 

2500 

6 

SO 

/ 

51 

36 

^o 

J  • 

1300 

6 

30 

/ 
/ 

35 

TO 

1300 

12 

30 

53 

23 

40 

3400 

10 

30 

/ 
/ 

54 

40 

2000 

30 

/ 

55 

40 

41 

1500 

10 

30 

/ 

56 

33 

41 

2600 

6 

80 

/ 

57 

42 

41 

1200 

10 

30 

/ 

53 

39 

42 

3500 

10 

SO 

5? 

39 

44 

3400 

6 

SO 

/ 

60 

42 

43 

1700 

12 

30 

/ 

61 

43 

44 

1400 

12 

30 

62 

45 

43 

750 

12 

SO 

63 

25 

46 

300 

6 

80 

/ 

64 

40 

46 

3000 

6 

SO 

/ 

65 

O  L 

■i  / 

1400 

16 

80 

/ 

66 

46 

47 

600 

6 

80 

/ 
/ 

67 

"Y  / 

48 

600 

16 

30 

68 

48 

49 

2500 

16 

80 

• 

6? 

49 

50 

1000 

12 

30 

/ 

70 

41 

50 

5000 

12 

80 

/ 

71 

50 

51 

3300 

16 

80 

/ 

7^ 

44 

51 

2000 

12 

80 

73 

48 

52 

2800 

16 

80 

/ 
/ 
/ 

74 

52 

53 

1100 

16 

80 

-»rr 

53 

54 

1600 

6 

30 

/ 

76 

49 

54 

1000 

30 

/ 

77 

53 

ir  tr 

JJ 

1100 

16 

30 

78 

55 

JO 

1800 

12 

30 

/ 

7? 

54 

56 

300 

12 

80 

/ 

SO 

5o 

J  / 

2500 

16 

30 

/ 

81 

58 

J  / 

1180 

24 

80 

32 

51 

cr.O 
J  u 

1600 

16 

30 

/ 

33 

57 

ero 
J  / 

3200 

20 

30 

/ 

84 

59 

60 

**♦  —**  f\.   ^ 

10 

80 

35 

58 

60 

3300 

o 

30 

q  l 

Jo 

rr  rr 

JJ 

61 

4100 

16 

30 

/ 

KJ   1 

9  /. 

4100 

16 

80 

/ 

38 

52 

2900 

6 

80 

/ 

89 

62 

A3 

1300 

12 

30 

/ 
/ 

90 

63 

64 

1800 

12 

80 

/ 
/ 

91 

/  4 

/  rr 

500 

12 

30 

/ 
/ 

700 

it 
o 

30 

i 

66 

3700 

6 

30 

/ 
/ 

94 

6o> 

61 

4900 

i.1 

30 

/ 

/ 

/  rr 

/.  "7 
J  / 

1200 

12 

30 

/ 

?o 

t  q 

1500 

12 

80 

/ 

,-j  -. 

£  "T 

•*▼  f\    r\    '^ 
\.s   \.y  v/  w 

/ 

Q  O 

/  ~7 

/  r\ 

a  / 

2000 

/ 
ij 

so 

/ 

no 

70 

.(.  o 

2200 

O 

30 

i 

00 

/  rr 

O  J 

70 

2  0  0  0 

c. 
o 

80 

/ 

53 


Table    13    (continued) 


NODE 

5 

1 

A-   » .1 

161 

n 

Am 

0 

302 

3 

13 

o  "Jf  cr 

/ 
/ 

A 

27 

/ 
t 
/ 

5 

0 

282 

/ 

6 

29 

/ 

7 

15 

120 

/ 

8 

27 

73 

9 

53 

/ 

10 

0 

55 

/ 

11 

10 

81 

/ 

12 

-650 

56 

/ 

13 

0 

55 

/ 

14 

0 

55 

/ 

15 

0 

55 

/ 

16 

0 

260 

/ 

17 

-1250 

55 

/ 

18 

91 

99 

/ 

1? 

63 

88 

/ 

20 

16 

96 

/ 

21 

7 

102 

/ 

9 

101 

/ 

23 

26 

106 

/ 

24 

34 

112 

/ 

25 

5 

120 

2  A 

21 

110 

/ 

27 

0 

84 

/ 

28 

0 

126 

/ 

29 

0 

210 

/ 

30 

34 

124 

/ 

31 

16 

127 

/ 

"SO 

20 

145 

/ 

j  3 

126 

69 

34 

13 

75 

/ 

35 

14 

62 

/ 

36 

30 

60 

/ 

37 

0 

277 

/ 

38 

0 

56 

39 

30 

70 

/ 

40 

15 

139 

41 

"to 

90 

42 

6 

**7  **"? 

/ 

43 

0 

47 

44 

42 

87 

/ 
t 

45 

0 

44 

/ 

46 

6 

137 

* 

3 

144 

43 

12 

130 

/ 

49 

— Y 

S 1 

/ 

50 

3 1 

56 

/ 

54 


Table    13    (continued) 
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v->  A.. 

53 

54 

55 

56 

57 

58 

5? 

60 

6.1 

Q<&. 

63 

64 

65 

66 

67 

68 

i  o 

70 


48 

4 

4 
'■\ 

A*. 

it  r\ 

>J ».') 

46 

41 

8 

26 
n  /. 

A„  l_> 

10 

r-\    * 

4 

? 

cr 

8 


..J 

13 


RESERVOIRS 


A« 

ir 


1  o 

37 


71'! 

w->    J.    A'.. 
\.)     J.     A— 

n  /.  n 

A--  I.)  AU 

250 
/->  — f  — . 

aC   /      / 


95 

101 

98 

91 
96 

8.1 

<■  ^ 

105 
60 
120 
67 
101 
116 
112 
111 
108 
116 
118 
109 
100 


SOURCE  PUMPS 
10  500  234 
14   200  360 

BOOSTER  PUMPS 
6   600  130 

CHECK  VALVES 


/ 

/ 
/ 
/ 


O       / 

CHECK 
END 

EH'D    OF 


600  264 

700 

OOcr 

A—  A..  »..' 

40  / 

J  0  (/   \j  ...  o 

400 

A.OJ 

40  / 

700  115 

800 

So 

0*0 

i::-  T 


LE 
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Much  of   this   data   was   subsequently   corrected   since   errors    in   the   data 
collection  process   had   not   then   been    identified.      However,    only    two 
"adjustments"   were  made   during   the  calibration  process,    one  being   the  use  of 
a  scaling    factor    for    the   user   demand   and    the   other   being    the   changing   of   the 
pipe  roughness   coefficients.      These   "adjustments"   were  made   to  get  agreement 
between    the   computer   model   and    the   ISO  calibration  data   set. 

The   1984   user   demand   estimated   in  section   II   totalled   only   4,000 
gallons   per  minute    or  5.76   million   gallons    per  day.      The   well    and   MDC 
records   indicated   that   during   the  calibration  period   the   water   supply   from 
the   wells   and   the  MDC   was   either   4,000   gpm  or   4,900   gpm  depending   on  whether 
well  C0   was  off  or  on.      Also  the   reservoir   level   records   showed   that   the 

c. 

Shaker  Glen  reservoirs  were  filling  slowly,  while  the  Rag  Rock  tank  had  a 
substantial  outflow,  roughly  700  gpm.   This  Rag  Rock  outflow  indicated  that 
the  user  demand  exceeded  the  well  and  MDC  supply  during  the  test  period. 
Therefore  a  scaling  factor  was  used  to  increase  the  demand  uniformly 
throughout  the  city.   First  the  demands  were  multiplied  by  a  scaling  factor 
of  1.5.   This  seemed  too  large  so  1 . 4  was  tried  and  found  to  be  in 
reasonable  agreement  with  the  calibration  data.   A  comparison  between  the 
supply  and  reservoir  flows  of  the  computer  model  and  the  calibration  data  is 
presented  in  Table  14. 

The  total  head  data  based  on  the  ISO  measurements  are  given  in  Table 
15.   They  show  that  when  the  hydrants  are  not  flowing  the  total  head  is 
roughly  250  feet  in  most  of  the  city  of  Woburn.   The  total  head  of  the 
computer  model  was  also  roughly  250  feet  from  the  first  trial,  with  C  =  80 
for  all  pipes,  until  the  final  roughness  values  were  determined.   The  first 
trial  produced  reasonable  total  heads,  but  it  also  showed  a  strong  outflow 
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Table  14.   Calibration  Flow  Comparison 
Supply  Computer  Model  Data 

Wells  A  ,  B,  Cp  (on) 

Well  D 
Well  E 
Well  F 
MDC 


Reservoir  Inflow 

Zion  Hill 
Whispering  Hill 
Rag  Rock 
Horn  Pond 


2208  gpm 

2208  gpm 

400 

400 

501 

454 

235 

262 

1817 

1604 

I 

=  5161 

I  =  4928 

Computer 

Model 

Data 

42  gpm 

100  gpm 

174 

None 

-459 

-625  (avg.) 

-123 

None 
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Table  15.   Total  Head  Based  on  1933  ISO  Data 


Location  Node/Pipe       P/Y   Elevation     H  (ft) 

(ft)    (ft) 


Cranes  Court  off  Main  P-23  170.8  77.5  248.3 

Main  and  Montvale  N-21  143.1  106.0  249.1 

Henshaw  and  Erie  P-50  216.9  57.0  273.9 

Montvale  and  Campbell  P-42  124.6  127.3  251.9 

Montvale  Road  and  Dale  P-58  173.1  78.2  251.3 

Dragon  Court  and  Derry  N-60  124.6  122.0  246.6 

Commerce  and  Mishawum  N-57  186.9  62.0  248.9 

Cedar  and  Washington  N-44  163.8  82.0  245.8 

Micro  and  Holton  P-45  189.2  74.0  263.2 

Commerce  and  Atlantic  N-59  175.4  75.2  250.6 

School  and  New  Boston  N-53  159.2  97.7  256.9 

New  Boston  and  Sixth  P-86  170.8  68.0  238.8 

New  Boston  §   end  N-61  170.8  71.0  241.8 

Main  and  Mountain  N-70  133.8  98.0  231.8 

Webster  St.  and  Webster  Ave.  N-67  120.0  119.0  239.0 

Middlesex  Canal  and  Main  P-87  131.5  115.0  246.5 

Main  and  Clinton  N-25  124.6  121.0  245.6 

Continental  Court  and  Salem  N-40  101.5  140.0  241.5 

Mishawum  and  Olympia  N-49  159.2  81.0  240.2 

Bedford  and  Houghton  P-1 3  138.4  110.2  248.6 

Bedford  and  Cambridge  N-7  126.9  128.0  254.9 

Totman  and  Lexington  P-9  177.7  53-0  230.7 

Russell  and  Silvermine  P-1  113.1  182.5  295.6 

Carol  and  Anna  N-1  138.4  159.6  293.0 

Lexington  and  Cambridge  N-4  242.3  71.0  313-3 

Frances  and  Main  N-24  120.0  123.5  243.6 

Cambridge  and  Country  Club   P-1 4  147.7  88.5  236.2 

Foster  and  Jones  P-94  147.7  98.0  245.7 
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from  the  Horn  Pond  reservoir  and  inflow  to  the  Rag  Rock  tank.   To  reduce 
this  transfer  flow  required  an  increased  pipe  roughness,  so  the  roughness  of 
the  main  pipes  was  changed  to  a  value  based  on  the  age  of  the  pipes.   This 
increase  brought  about  the  desired  outflow  from  the  Rag  Rock  tank,  but  it 
lowered  the  total  head  to  roughly  240  feet  and  lowered  the  head  in  northwest 
Woburn  well  below  the  calibration  values.   Thus  the  roughness  along  Main 
Street  north  of  City  Hill  was  returned  to  C  =  80.   This  change  produced  a 
computer  model  whose  total  head  and  flow  distributions  were  close  to  the 
calibration  data. 

The  computer  model  was  then  compared  with  the  measurements  taken  with 
the  hydrants  flowing.   This  data  is  presented  later,  in  the  section  on  the 
error  analysis  of  this  calibration  process  (Table  17b).   Generally  these 
hydrant  tests  indicated  that  the  pipe  roughness  was  too  high,  and  so  the  C 
values  were  raised  as  needed.   The  resulting  "calibrated"  flow  model  for 
the  conditions  of  the  ISO  measurement  days  is  presented  in  Table  16. 

A  flow  map  for  the  calibrated  flow  model  is  presented  in  Figure  16. 
The  noteworthy  features  of  this  flow  map  are  the  south  to  north  general  flow 
pattern  due  to  the  southerly  location  of  the  Woburn  wells,  the  east  to  west 
general  flow  pattern  due  to  the  total  head  of  the  MDC  connection  and  the  low 
total  head  of  the  suction  side  of  the  Shaker  Glen  booster  pump,  and  the 
boundary  between  the  Woburn  well  water  and  the  MDC  water. 

The  north-south  flow  is  well  defined  and  easily  seen  due  to  the 
difference  in  water  elevation  between  the  Horn  Pond  and  Rag  Rock  reservoirs. 
The  12  foot  head  drop  gives  a  clear  indication  of  the  north-south  energy 
gradient  which  continues  along  Main  Street.   Even  when  Well  C?  is  turned 

off,  the  northerly  flow  continues  strong  since  the  Horn  Pond  reservoir  has 
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Table    16 .      Calibrated    Model   - 


PIPE  DATA: 

Compu 

ter  Output 

PIPE 

NODES 

HEAD 

HLOSS 

NO. 

FROM 

TO 

LENGTH 

DIAM 

COEF 

FLOW  RATE 

VELOCITY 

LOSS 

/1000 

1 

4 

1 

4200. 

6.0 

40.0 

64.04 

.73 

17.02 

4.05 

2 

3 

1 

5200. 

10.0 

40.0 

197.76 

.81 

14.13 

2.72 

3 

4 

3 

3700. 

12.0 

90.0 

366.22 

1.04 

2.89 

.78 

»   4 

3 

2 

700. 

10.0 

90.0 

42.46 

.17 

.02 

.04 

*   5 

4 

5 

4600. 

12.0 

90.0 

174.61 

.50 

.91 

.20 

6 

6 

4 

400. 

12.0 

90.0 

753.27 

2.14 

1.19 

2.97 

*   7 

7 

6 

4800. 

6.0 

80.0 

80.89 

.92 

8.31 

1.73 

8 

8 

7 

2800. 

10.0 

80.0 

233.49 

.95 

2.87 

1.02 

9 

9 

6 

3000. 

6.0 

40.0 

70.00 

•  79 

14.34 

4.78 

10 

9 

8 

4500. 

12.0 

80.0 

307.86 

.87 

3.16 

.70 

11 

10 

9 

800. 

12.0 

50.0 

501.15 

1.42 

3.31 

4.14 

*   12 

11 

9 

2600. 

6.0 

40.0 

18.11 

.21 

1.02 

.39 

13 

24 

8 

7000. 

10.0 

90.0 

90.83 

.37 

1.00 

.14 

14 

12 

6 

7200. 

12.0 

90.0 

789.97 

2.24 

23.32 

3.24 

15 

12 

13 

5700. 

16.0 

90.0 

716.53 

1.14 

3.80 

.67 

16 

13 

11 

3700. 

16.0 

90.0 

951.66 

1.52 

4.17 

1.13 

17 

14 

13 

100. 

10.0 

90.0 

235.13 

.96 

.08 

.83 

*     18 

-15 

12 

800. 

14.0 

80.0 

1106.10 

2.31 

2.84 

3.55 

19 

16 

15 

500. 

16.0 

40.0 

122.67 

.20 

.06 

.11 

*  20 

17 

•  15 

200. 

16.0 

40.0 

983.43 

1.57 

1.07 

5.37 

21 

17 

18 

6000. 

16.0 

40.0 

578.68 

.92 

12.08 

2.01 

22 

11 

18 

1300. 

16.0 

40.0 

145.07 

.23 

.20 

.16 

23 

17 

19 

4700. 

16.0 

40.0 

645.69 

1.03 

11.59 

2.47 

24 

19 

20 

900. 

16.0 

40.0 

394.98 

•  63 

.89 

.99 

25 

20 

21 

2000. 

16.0 

40.0 

238.53 

•  38 

.78 

•  39 

26 

21 

22 

300. 

16.0 

40.0 

123-50 

.20 

.03 

.12 

27 

18 

22 

1300. 

16.0 

40.0 

383.10 

-.61 

1.22 

.94 

28 

22 

23 

900. 

16.0 

80.0 

821.18 

1.31 

.96 

1.07 

29 

23 

24 

1400. 

16.0 

80.0 

585.32 

.93 

.80 

.57 

30 

24 

25 

1900. 

16.0 

80.0 

238.29 

.38 

.20 

.11 

31 

25 

26 

700. 

16.0 

80.0 

188.72 

•  30 

.05 

.07 

32 

11 

27 

1700. 

16.0 

90.0 

686.28 

1.10 

1.05 

.61 

33 

.18 

27 

1100. 

10.0 

40.0 

99.85 

.41 

.84 

.77 

i  34 

27 

22 

2100. 

16.0 

90.0 

352.38 

.56 

.38 

.18 

35 

27 

28 

1750. 

20.0 

90.0 

433.75 

.44 

.16 

.09 

36 

29 

28 

1400. 

24.0 

90.0 

458.82 

.33 

.06 

.04 

37 

28 

26 

3800. 

16.0 

120.0 

892.58 

1.42 

2.23 

.59 

38 

19 

30 

1500. 

8.0 

40.0 

63. 11 

.40 

1.46 

.97 

39 

20 

30 

1700. 

10.0 

50.0 

79.44 

•  32 

.56 

.33 

40 

21 

31 

1100. 

12.0 

120.0 

63-23 

.18 

.02 

.02 

41 

30 

31 

3400. 

8.0 

60.0 

22.70 

.14 

.23 

.07 

t     42 

32 

31 

1900. 

12.0 

60.0 

73.66 

.21 

.16 

.08 

43 

33 

32 

2200. 

12.0 

60.0 

224.65 

.64 

1.47 

.67 

44 

34 

33 

1000. 

12.0 

60.0 

546.93 

1.55 

3.47 

3.47 

45 

34 

30 

4500. 

10.0 

50.0 

150.35 

.61 

4.87 

1.08 

46 

35 

34 

1500. 

12.0 

60.0 

764.48 

2.17 

9.69 

6.46 

47 

36 

35 

600. 

12.0 

60.0 

1563.76 

4.44 

14.58 

24.30 

48 

37 

36 

100. 

12.0 

120.0 

1816.89 

5.15 

.89 

8.89 
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Table 

16  (cont 

inued  ) 

PIPE  DATA: 

PIPE 

NODES 

HEAD 

HLOSS 

NO. 

FROM 

TO 

LENGTH 

DIAM 

COEF 

FLOW  RATE 

VELOCITY 

LOSS 

/1000 

»   4g 

38 

35 

3300. 

6.0 

40.0 

46.37 

.53 

7.35 

2.23 

50 

36 

38 

2500. 

6.0 

40.0 

53-37 

.61 

7.23 

2.89 

51 

36 

39 

1800. 

6.0 

50.0 

126.97 

1.44 

17.14 

9.52 

52 

35 

39 

1300. 

12.0 

120.0 

805.05 

2.28 

2.56 

1.97 

53 

23 

40 

3400. 

10.0 

120.0 

104.25 

•  43 

.37 

.11 

54 

32 

40 

2000. 

6.0 

80.0 

51.59 

.59 

1.50 

.75 

*  55 

41 

40 

1500. 

10.0 

120.0 

10.96 

.04 

.00 

.00 

56 

33 

41 

2600. 

6.0 

80.0 

64.68 

.73 

2.97 

1.14 

57 

42 

41 

1200. 

10.0 

80.0 

243.41 

.99 

1.33 

1.11 

58 

39 

42 

3500. 

10.0 

110.0 

623.28 

2.55 

12.25 

3.50 

59 

39 

44 

3400. 

6.0 

50.0 

80.53 

.91 

13.93 

4.10 

60 

42 

43 

1700. 

12.0 

100.0 

335.08 

.95 

.93 

.54 

61 

43 

44 

1400. 

12.0 

100.0 

335.08 

.95 

.76 

.54 

62 

45 

43 

750. 

12.0 

100.0 

.00 

.00 

.00 

.00 

63 

25 

46 

800. 

6.0 

50.0 

22.97 

.26 

.32 

.40 

64 

40 

46 

3000. 

6.0 

80.0 

32.41 

•  37 

.95 

.32 

65 

26 

47 

1400. 

16.0 

90.0 

368.50 

.59 

.27 

.19 

*  66 

.47 

46 

600. 

6.0 

90.0 

.62 

.01 

.00 

.00 

67 

47 

48 

600. 

16.0 

90.0 

342.68 

.55 

.10 

.17 

«  68 

49 

48 

2500. 

16.0 

90.0 

23.17 

.04 

.00 

.00 

«  69 

50 

.  49 

1000. 

16.0 

90.0 

66.45 

.11 

.01 

.01 

70 

41 

50 

5000. 

12.0 

110.0 

220.12 

.62 

1.05 

.21 

71 

50 

51 

3800. 

16.0 

90.0 

114.46 

.18 

.08 

.02 

72 

44 

51 

2000. 

12.0 

100.0 

278.41 

.79 

.77 

.39 

73 

48 

52 

2800. 

16.0 

110.0 

294.45 

.47 

.25 

.09 

74 

52 

53 

1100. 

16.0 

110.0 

208.26 

.33 

.05 

.05 

*  75 

54 

53 

1600. 

6.0 

50.0 

3.16 

-.04 

.02 

.01 

76 

49 

54 

1000. 

6.0 

40.0 

15.29 

.17 

.29 

.29 

77 

53 

55 

1100. 

16.0 

110.0 

108.09 

.17 

.02 

.01 

*  78 

56 

55 

1800. 

12.0 

110.0 

57.80 

-.16 

.03 

.02 

*  79 

56 

54 

800. 

12.0 

110.0 

6.07 

.02 

.00 

.00 

»  80 

57 

56 

2500. 

16.0 

120.0 

128.27 

.20 

.04 

'   .02 

81 

58 

57 

1100. 

24.0 

110.0 

264.98 

.19 

.01 

.01 

82 

51 

58 

1600. 

16.0 

110.0 

311.67 

.50 

.16 

.10 

83 

57 

59 

3200. 

20.0 

110.0 

72.30 

.07 

.01 

.00 

84 

59 

60 

3700. 

10.0 

120.0 

34.50 

.14 

.05 

.01 

85 

58 

60 

3300. 

6.0 

50.0 

4.70 

.05 

.07 

.02 

86 

55 

61 

4100. 

16.0 

110.0 

83-30 

.13 

.03 

.01 

87 

26 

62 

4100. 

16.0 

80.0 

513.99 

.82 

1.84 

.45 

88 

52 

62 

2900. 

6.0 

100.0 

46.99 

.53 

1.21 

.42 

89 

62 

63 

1300. 

12.0 

80.0 

489.58 

1.39 

2.16 

1.66 

90 

63 

64 

1800. 

12.0 

120.0 

402.08 

1.14 

.98 

.54 

91 

64 

65 

500. 

12.0 

120.0 

422.10 

1.20 

.30 

.60 

»  92 

66 

64 

700. 

6.0 

90.0 

49.42 

.56 

.39 

.56 

93 

53 

66 

3700. 

6.0 

90.0 

69.73 

.79 

3-91 

1.06 

i  94 

61 

66 

4900. 

6.0 

90.0 

59.50 

.68 

3.86 

.79 

95 

65 

67 

1200. 

12.0 

120.0 

297.68 

.84 

.37 

.31 

96 

67 

68 

1500. 

12.0 

80.0 

168.70 

.48 

-35 

.23 

97 

63 

68 

3000. 

6.0 

80.0 

48.30 

.55 

2.00 

.67 

98 

67 

69 

2000. 

6.0 

120.0 

70.18 

.80 

1.26 

.63 

»  99 

69 

70 

2200. 

6.0 

120.0 

22.58 

.26 

.17 

.08 

100 

65 

70 

2000. 

6.0 

120.0 

85.22 

.97 

1.80 

.90 

61 


?able    lo    (continued) 


NODE 

DEMAND 

HGL 

NO. 

(CFS) 

(GPM) 

ELEV 

HEAD 

PRESSURE 

ELEV 

1 

.583 

261.80 

161.00 

135.89 

58.89 

296.89 

3 

.281 

126.00 

235.00 

76.02 

32.94 

311.02 

4 

•  331 

148.40 

64.00 

249.91 

108.29 

313-91 

6 

.418 

187.60 

62.00 

170.79 

74.01 

232.79 

7 

.340 

152.60 

120.00 

121.10 

52.47 

241.10 

8 

.368 

165.20 

73.00 

170.96 

74.08 

243.96 

9 

•  315 

141.40 

53-00 

194.13 

84.12 

247.13 

11 

.228 

102.20 

81.00 

167.14 

72.43 

248.14 

12 

-.892 

-400.40 

56.00 

200.11 

86.71 

256.11 

13 

0.000 

0.00 

47.00 

205.31 

88.97 

252.31 

15 

0.000 

0.00 

60.00 

198.94 

86.21 

258.94 

17 

-4.919 

-2207.80 

47.00 

213.02 

92.31 

260.02 

18 

.537 

240.80 

99.00 

148.94 

64.54 

247.94 

19 

.418 

187.60 

88.00 

160.43 

69.52 

248.43 

20 

.172 

77.00 

96.00 

151.54 

65.67 

247.54 

21 

.115 

51.80 

102.00 

144.76 

62.73 

246.76 

22 

.084 

37.80 

101.00 

145.72 

63-15 

246.72 

23 

.293 

131.60 

106.00 

139.76 

60.56 

245.76 

24 

.571 

256.20 

112.00 

132.97 

57.62 

244.97 

25 

.059 

26.60 

120.00 

124.76 

54.06 

244.76 

26 

.443 

198.80 

110.00 

134.71 

58.38 

244.71 

27 

0.000 

0.00 

84.00 

163.10 

70.68 

247.10 

28 

0.000 

0.00 

126.00 

120.94 

52.41 

246.94 

30 

.602 

270.20 

124.00 

122.97 

53.29 

246.97 

31 

.356 

159.60 

127.00 

119.74 

51.89 

246.74 

32 

.221 

99.40 

145.00 

101.90 

44.16 

246.90 

33 

.574 

257.60 

69.00 

179.37 

77.73 

248.37 

34 

.150 

67.20 

75.00 

176.84 

76.63 

251.84 

35 

.090 

40.60 

62.00 

199.53 

86.46 

261.53 

36 

.162 

72.80 

60.00 

216.11 

93.65 

276.11 

38 

.016 

7.00 

56.00 

212.88 

92.25 

268.88 

39 

.508 

228.20 

70.00 

188.97 

81.89 

258.97 

40 

.299 

134.40 

139.00 

106.39 

46.10 

245.39 

41 

.172 

77.00 

90.00 

155.40 

67.34 

245.40 

42 

.100 

44.80 

77.00 

169.72 

73.55 

246.72 

43 

0.000 

0.00 

47.00 

198.80 

86.15 

245.80 

44 

.306 

137.20 

87.00 

158.04 

68.48 

245.04 

45 

0.000 

0.00 

44.00 

201.80 

87.45 

245.80 

46 

.125 

56.00 

137.00 

107.44 

46.56 

244.44 

47 

.056 

25.20 

144.00 

100.44 

43.52 

244.44 

48 

.159 

71.40 

130.00 

114.34 

49.55 

244.34 

49 

.062 

28.00 

81.00 

163.34 

70.78 

244.34 

50 

.087 

39.20 

56.00 

188.35 

81.62 

244.35 

51 

.  181 

81.20 

95.00 

149.26 

64.68 

244.26 

52 

.087 

39.20 

101.00 

143.09 

62.01 

244.09 

53 

.075 

33-60 

98.00 

146.04 

63.28 

244.04 

54 

.041 

18.20 

91.00 

153.06 

66.32 

244.06 

55 

.184 

82.60 

96.00 

148.02 

64.14 

244.02 

56 

.143 

64.40 

81.00 

163-06 

70.66 

244.06 

57 

.143 

64.40 

62.00 

182.10 

78.91 

244. 10 

58 

.094 

42.00 

105.00 

139.11 

60.28 

244.11 

59 

.084 

37.80 

60.00 

184.09 

79.77 

244.09 

60 

.087 

39.20 

120.00 

124.04 

53.75 

244.04 

61 

.053 

23-80 

67.00 

176.99 

76.70 

243.99 

62 

.159 

71.40 

101.00 

141.88 

61.48 

242.88 

63 

.087 

39.20 

1 1 6 . 00 

124.72 

54.04 

240.72 

64 

.066 

29.40 

112.00 

127.74 

55.35 

239.74 

65 

.087 

39.20 

111.00 

128.44 

55-66 

239.44 

66 

.178 

79.80 

108.00 

132.13 

57.26 

240.13 

67 

•  131 

58.80 

1 1 6 . 00 

123.07 

53-33 

239-07 

68 

.483 

217.00 

118.00 

120.72 

52.31 

238.72 

69 

.106 

47.60 

109.00 

128.81 

55.82 

237.81 

70 

.240 

107.80 

100.00 

137.64 

59.65 

237.64 

10 

-1.117 

-501.15 

43.00 

207.44 

89.89 

250.44 

14 

-.524 

-235.13 

46.00 

206.39 

89.44 

252.39 

2 

.095 

42.46 

302.00 

9.00 

3-90 

311.00 

5 

.389 

174.61 

282.00 

31.00 

13.43 

313-00 

16 

-.273 

-122.67 

239.00 

20.00 

8.67 

259.00 

29 

-1.022 

-458.82 

210.00 

37.00 

16.03 

247.00 

37 

-4.048 

-1816.89 

277.00 

0.00 

0.00 

277.00 
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almost  the  same  total  head  as  does  well  C2,  and  so  the  reservoir  outflow 
easily  replaces  the  water  from  well  C  . 

The  44  foot-head  drop  between  the  MDC  connection  and  the  Shaker  Glen 
booster  pump  at  Russell  and  Cambridge  Streets  causes  a  westerly  flow  along 
Bedford  Road  and  even  brings  the  MDC  water  southwest  along  Mishawum  Road  as 
far  as  Forest  Park  Road.   There  it  is  able  to  return  to  its  westerly  path 
into  northwest  Woburn.   When  this  booster  pump  is  turned  off,  there  is  a 
counterflow  from  west  to  east  due  to  well  E,  but  when  it  is  on,  the  booster 
pump  easily  dominates  well  E.   The  boundary  between  the  Woburn  well  water 
and  the  MDC  water  is  identified  by  the  heavy  line  running  roughly  north- 
south  just  east  of  Main  Street.   Mixing  of  water  from  the  two  sources  occurs 
at  nodes  #30,  40,  and  48.   These  mixtures  then  move  north  along  pipes  which 
carry  the  mixed  water  away  from  the  mixing  pipe  junctions.   These  pipes, 
numbered  41,  64,  and  73  respectively,  form  a  part  of  the  boundary  shown  in 
Figure  16.   Even  stronger  mixing  occurs  in  northwest  Woburn.   Since  pipe  #73 
carries  mixed  water  to  both  Alfred  Street  and  New  Boston  Street  all  the 
water  in  northwest  Woburn  is  a  mixture  of  water  from  the  two  sources. 

The  error  level  of  this  ISO  calibration  process  is  obtained  by 
comparing  the  total  head  distribution  of  the  ISO  measurements  with  that  of 
the  computer  model  after  all  the  "adjustments"  to  peaking  factor  and  pipe 
roughness  were  made.   Table  17a  presents  this  total  head  comparison  for  the 
flow  distribution  when  no  hydrants  are  flowing.   The  average  head  difference 
of  the  28  calibration  points  was  0.9  feet  or  0.4  psi,  while  the  root  mean 
square  head  difference  was  6.15  feet  or  2.7  psi.   Further  adjustment  seems 
unwarranted.   Table  17b  presents  a  comparison  of  the  pressure  measured  at 
the  hydrant  near  the  different  flowing  hydrants  and  the  pressure  predicted 
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Table  17a.   Calibration  Head  Error 


Node  // 

1 
H 

7 
21 
24 

25 

40 

44 

49 

53 

57 

59 

60 

61 

67 

70 


Data  H(ft) 

Model  H(ft) 

(HD"HM) 

298.0 

296.9 

2.1 

313.3 

313.9 

-0.6 

254.9 

241  .1 

13.8 

249.1 

246.8 

2.3 

243.5 

245.0 

-1  .5 

245.6 

244.8 

0.8 

241  .5 

245.4 

-3.9 

250.8 

245.0 

5.8 

240.2 

244.3 

-4.1 

256.9 

244.0 

12.9 

248.9 

244.1 

4.8 

250.6 

244.1 

6.5 

246.6 

244.0 

2.6 

241  .8 

244.0 

-2.2 

239.0 

239.1 

-0.1 

231  .8 

237.6 

-5.8 

Pipe  # 


H  (ft) 


H(ft) 


1 

295.6 

300.4 

9 

230.7 

241  .6 

13 

248.6 

244.5 

14 

236.2 

244.4 

23 

248.3 

254.2 

42 

251  .9 

246.8 

45 

263.2 

250.6 

50 

273.9 

274.0 

58 

251.3 

253-0 

87 

246.5 

243.8 

86 

238.8 

244.0 

94 

245.7 

242.0 

N    =   23 

Avg. 

=    0.89 

ft 

S 

=   6.15 

ft 

-4.8 

-10.9 

4.1 

-8.2 

-5.9 

5.1 

12.6 

-0.1 

-1.7 

2.7 

-5.2 

3.7 
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Table  17b.   Pressure  Error  for  Hydrant  Flows 
Node  No.  Data  (psi) 


7  18 

21  59 

25  50 

31  51 

40  36 

44  58 

49  67 

53  55 

57  73 

61  36 

67  26 

70  25 

N  =  12  Avg.  =  1 .9  psi 

S  =  15.4  psi 


Computer   Model 

AP 

(psi) 

-17.4 

35 

59.1 

0 

48.3 

2 

47.8 

3 

36.3 

0 

57.2 

1 

62.6 

4 

57.7 

-3 

75.3 

-2 

69.9 

-34 

23.7 

2 

10.2 

15 
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by  the  computer  model  at  the  nearby  pipe  junctions.   The  average  pressure 
difference  is  1.9  psi  while  the  root  mean  square  difference  is  15.4  psi. 
This  large  root  mean  square  difference  is  due  to  three  cases  where  the 
agreement  between  the  model  and  the  data  is  quite  poor.   At  node  #7,  Bedford 
Road  and  Cambridge  Street,  the  model  predicts  a  much  lower  pressure  than 
that  measured.   The  model  prediction  assumed  the  Shaker  Glen  booster  pump 
was  on,  while  the  status  of  the  pump  when  the  measurement  was  made  is 
unknown.   Since  the  local  C  values  are  already  high,  the  only  reasonable 
cause  of  this  error  is  that  the  booster  pump  "must"  have  been  temporarily 
off  when  the  measurement  was  taken.   At  node  #61 ,  New  Boston  Street  and 
Merrimac  Street  (see  Figure  16),  the  model  predicts  a  much  higher  pressure 
than  that  measured.   The  reason  is  that  the  site  of  the  measurement  was  at 
the  end  of  Mew  Boston  Street,  1800  feet  of  six  inch  pipe  away  from  the  16 
inch  pipe  at  New  Boston  and  Merrimac.   There  is  a  large  pressure  drop 
required  to  get  the  1290  gpm  fire  flow  through  this  "extra"  pipe.   The 
situation  at  node  #70,  Main  and  Mountain  Streets,  is  not  so  easily 
explained.   The  model  prediction  is  much  lower  than  the  measured  pressure. 
The  local  values  of  C  are  already  as  high  or  even  higher  than  is  reasonable. 
The  hydrant  flow  was  measured  and  so  the  demand  is  probably  correct.   The 
pipe  skeletonization  was  checked  and  also  compared  with  that  of  Waldorf  and 
Cleary  and  no  error  was  found.   With  the  data  available  the  cause  of  this 
error  cannot  be  detected.   Thus  two  of  the  three  problem  hydrant  tests  can 
be  explained  but  not  the  third.   Neglecting  those  two  the  average  pressure 
difference  is  2.2  psi  and  the  root  mean  square  difference  is  5.0  psi,  with 
the  latter  value  due  mainly  to  node  #70. 
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IV.   Computer  Model  Validation 

The  validity  of  the  computer  model  was  checked  by  field  tests  of  the 
pressure  and  fluoride  concentration  distributions  in  Woburn.   The  field 
tests  were  a  joint  effort  of  the  Woburn  water  Department,  the  DEQE  Water 
Supply  Division,  and  the  author.   Two  different  model  predictions  were 
checked,  the  first  being  the  pressure  distribution  near  Cedar  and  Washington 
Streets  when  the  hydrant  at  that  intersection  was  flowing  and  the  second 
being  the  pressure  and  fluoride  distributions  throughout  Woburn  when  no 
hydrants  were  flowing.   These  two  field  tests  occurred  on  the  afternoon  of 
June  4,  1985,  which  was  arbitrarily  chosen  as  a  typical  day.   A  secondary 
test  that  was  not  compared  with  a  model  prediction  was  the  early  morning  (on 
June  4,  1981  between  6  and  7  AM)  measurement  of  the  fluoride  concentration 
along  the  pipes  that  connect  the  MDC  fluoride  source  with  the  Rag  Rock  tank. 

The  computer  model  prediction  is  based  on  the  assumptions  that  November 
2  and  3,  1983  and  June  4,  1985  were  all  typical  days  of  water  demand, 
supply,  and  distribution  for  the  City  of  Woburn,  and  that  the  hydrant  would 
be  opened  the  same  way  during  the  validation  test  as  it  was  during  the  ISO 
tests.   If  major  differences  in  the  water  system  operation  during  the  two 
tests  were  noted,  then  the  computer  model  would  be  changed  to  account  for 
those  differences.   However  the  computer  model,  calibrated  for  the  typical 
days  in  1983,  would  not  be  "adjusted"  again;  the  pipe  roughness  and  user 
demand  distribution  would  not  be  changed.   Only  operational  changes,  like 
turning  pumps  on  or  off,  were  allowed  during  this  validation  calculation 
with  the  computer  model.   In  fact  a  validation  calculation  was  made  since 
major  differences  were  noted. 

Considerable  care  was  taken  to  insure  the  accuracy  of  the  field  tests. 
The  afternoon  period  from  1  PM  to  5  PM  was  chosen  as  the  period  during  the 
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day  with  the  most  steady  flow  conditions,  and  the  data  was  collected  as 
quickly  as  possible.   Both  these  considerations  were  made  to  assure  the 
validity  of  the  model-data  comparison.   The  pressure  gauges  used  in  the  test 
were  calibrated  by  comparison  with  the  Woburn  Pump  Station  pressure  gauge  on 
the  Well  C  pipe  and  with  the  Horn  Pond  reservoir  water  level.   The 

calibration  data  is  presented  in  Table  18.   With  well  C?  on  during  the 

calibration,  the  reference  to  the  Horn  Pond  reservoir  elevation  was  poor 
since  there  is  an  unknown  head  loss  between  that  gauge  and  the  reservoir. 
Only  a  relative  comparison  among  the  pressure  gauges  could  be  made.   The 
root  mean  square  difference  between  the  individual  gauges  and  their  mean  was 
only  0.6  psi,  indicating  that  they  closely  agree  with  one  another.   The 
fluoride  samples  were  carefully  collected  in  specially  prepared  bottles 
supplied  by  the  Lawrence  Experimental  Station  of  DEQE.   The  samples  were 
analyzed  at  both  the  University  of  Massachusetts,  Amherst  and  the  Lawrence 
Laboratory.   Further,  at  the  University  two  different  laboratory  techniques 
were  used.   Thus  the  fluoride  concentrations  could  be  cross-checked  two  ways 
to  guarantee  their  accuracy. 

When  the  hydrant  test  was  run,  there  were  enough  major  water  system 
differences  that  a  validation  calculation  was  needed.   Table  19a  shows  the 
typical  day  prediction,  the  validation  result  and  the  actual  measurements. 
The  three  system  changes,  279  feet  of  head  instead  of  277  feet  at  the  MDC 
connection,  two  booster  pumps  on  instead  of  one  at  Shaker  Glen,  and  only 
1,680  gpm  instead  of  2,175  gpm  as  the  demand  at  the  node  corresponding  to 
the  flowing  hydrant,  caused  a  significant  difference  between  the  typical  day 
prediction  and  the  validation  result. 
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Table  18.   Pressure  Meter  Calibration 


Reading  at  "0"    Read  at  "Test" 

0.0  92.0-0.0  =  92.0 

0.0  93.5-0.0  =  93.5 

1  .0  93.8-1 .0  =  92.8 

0.0  92.5-0.0  =  92.5 

1  .0  94.0-1 .0  =  93.0 

1.5  93.5-1 .5  =  92.0 


N  =  6  Avg.   =  92.6  psi 

S   =  0.6  psi 


Gauge  No. 

Gauge  Type 

1 

March  Inst.  Co 

(100  +  2  psi) 

2 

U.S.  Gauge 

(100  +  1  psi) 

3 

Ashcrof t 

(100  +  1  psi) 

4 

U.S.  Gauge 

(100  +  1  psi) 

5 

Ashcrof t 

(160  +  1  psi) 

6 

U.S.  Gauge 

(160  +  1  psi) 
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The  validation  result  was  roughly  three  psi  below  the  measured  values 
when  the  hydrant  at  Cedar  and  Washington  was  flowing,  but  Table  20  also 
shows  roughly  that  same  difference  for  the  no-flow  data.   Thus  a  comparison 
in  Table  19b  of  the  pressure  drops  that  occurred  when  the  hydrant  was  opened 
does  not  show  this  systematic  prediction  error.   The  pressure  drop 
calculation  shows  a  mean  error  of  0.2  psi  and  a  root  mean  square  error  of 
1 .8  psi. 

The  results  of  the  pressure  tests  with  no  flowing  hydrants  are 
presented  in  Table  20.   The  validation  result  in  this  case  identified  fairly 
important  hydrant  elevation  errors,  which  were  corrected.   The  elevations  of 
the  test  hydrants  were  used,  instead  of  the  node  elevations,  in  calculating 
the  validation  hydrant  pressures.   The  resulting  pressures  have  a  mean  error 
of  0.06  psi  and  a  root  mean  square  error  of  3.3  psi. 

A  typical  areal  distribution  of  the  fluoride  concentration  in  Woburn's 
water  was  obtained  from  57  water  samples  taken  from  53  different  sites 
during  the  afternoon  of  June  H,    1985.   In  a  separate  test,  six  samples  were 
taken  between  6  and  7  AM  that  same  day  to  check  the  effect  of  the  daily 
demand  cycle.   The  water  at  the  sample  sites  was  run  for  roughly  three 
minutes  before  the  sample  was  taken  in  order  to  get  samples  from  the  water 
mains.   The  samples  were  collected  in  acid-washed  polyethylene  bottles 
provided  by  the  Lawrence  Experiment  Station  of  DEQE.   Each  sample  was  tested 
by  the  electrode  method,  both  at  the  University  of  Massachusetts 
Environmental  Engineering  Laboratory  and  at  the  Lawrence  Laboratory. 
Further,  roughly  half  of  the  samples  were  also  analyzed  by  the  SPADNS 
colorimetric  method  at  the  University.   The  University  analysis  was  done  by 


71 


Table  19a.  Validation  Test  with  Hydrant  Flowing  at  Cedar  and  Washington 


Location 


Typical  Day  Validation 
Node  No.    Data   Prediction   Calculation 


Cedar  and  Washington    44 
Olympia  and  Mishawum    49 


Salem  and  Pine         42 


Washington  and  Pine     39 


(psi)   (psi) 


64.0 


70.0 


72.0 


57.2 


65.5 


67.0 


83.0    79.3 


(psi) 

61.3 
67.0 
69.4 

80.6 


Validation 
Calculation- 
Data 
(psi) 

-2.7 
-3.0 

-3.1 
-2.4 


Table  19b.   Validation  Test:   Pressure  Drop  Analysis 


Location 


Node  No. 


Cedar  and  Washington   44 
Olympia  and  Mishawum   49 


Salem  and  Pine 


42 


Washington  and  Pine    39 


Pressure  Pressure  Drop 
Drop  (Data)  (Calculation) 
(psi)       (psi) 


Diffrence  in 
Pressure  Drop 
(psi) 


5.0 

7.2 

2.2 

6.0 

3.8 

-2.2 

4.5 

4.2 

-0.3 

2.0 

1  .7 

Avg. 
S 

-0.3 
=  -0.15   psi 
1 .8   psi 
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Table  20.   Validation  Test  with  Hydrants  not  Flowing 


Location 


Cedar  and  Washington 


Olympia  and  Mishawum 


Salem  and  Pine 


Washington  and  Pine 


Typical  Day  Validation 
Node  No.   Data   Prediction   Calculation 
(psi)   (psi)        (psi) 


44      69.0    68.5 


49      76.0    70.8 


42      76.5    73.5 


39      85.0    81.9 


68.5 
70.8 
73-6 
82.3 


Validation 

Calculation-Data 

(psi) 

-0.5 
-5.2 

-2.9 


-2.7 


and  Cambridge 
and  Cambridge 


Russell 

(High) 
Russell 

(Low) 
Bedford  and  Cambridge 
Arlington  and  Pleasant 
Main  and  Clinton 
Green  and  Garfield 
Montvale  and  Bow 
Montvale  and  Washington 
Salem  and  Wildwood 
New  Boston  and 

New  Industrial 
Mishawum  and  Commerce 
Elm  and  West 

N  -  16 


110.0   108.3 


108.7 


-1.3 


6 

63.0 

74.0 

71  .4 

8.4 

7 

58.0 

52.5 

52.8 

-5.2 

18 

63.0 

64.5 

64.5 

1.5 

25 

52.0 

54.0 

54.1 

2.1 

30 

47.0 

53.3 

47.2 

0.2 

32 

38.0 

44.2 

37.7 

-0.3 

36 

95.0 

93.6 

94.5 

-0.5 

41 

65.0 

67.3 

67.4 

2.4 

55 

61  .0 

64.1 

64.2 

3.1 

57 

71  .0 

78.9 

72.0 

1  .0 

63 

55.0 

54.0 

54.0 

Avg.  = 
S  = 

-0.06 
3-3 

-1  .0 
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Candace  Balmer,  a  master's  degree  student.   A  detailed  description  of  her 
analysis  procedure  is  included  in  Appendix  V. 

The  results  of  the  fluoride  tests  are  given  in  Table  21  and  are  shown 
in  Figures  17  and  18.   The  afternoon  tests  show  three  zones,  a  Woburn  water 
zone  in  southwest  Woburn  with  concentrations  less  than  0.10  parts  per 
million  (ppm),  a  MDC  water  zone  with  concentrations  above  0.90  ppm,  and  a 
mixing  zone  between  the  other  two  with  intermediate  concentrations.   It  is 
noteworthy  that  Main  Street  north  of  Green  and  also  the  Rag  Rock  tank  are  in 
this  mixed  water  zone.   The  morning  test  shows  roughly  the  same  east-west 
distribution  along  the  MDC  connection  -  Rag  Rock  tank  line.   The  mixed 
water  zone  extends  between  North  Warren  Street  and  Montvale  Avenue  near  Main 
Street. 

The  error  level  indicated  by  the  difference  between  the  Lawrence  and 
University  electrode  test  results  is  very  small  and  is  considered 
insignificant.   The  SPADNS  results  were  in  general  agreement  but  of  much 
lower  accuracy  and  so  were  neglected. 

The  computer  model's  validation  results  can  also  be  used  to  analyze  the 
fluoride  concentration  and  thus  get  a  "validation  calculation"  of  the 
fluoride  concentration.   Assuming  a  concentration  of  0.07  ppm  for  the  Woburn 
water  and  1 .07  ppm  for  the  MDC  water  gives  the  areal  distribution  presented 
in  Figure  19.   All  mixing  is  assumed  to  occur  at  nodes  where  the  mass 
conservation  principle  is  applied  to  calculate  the  concentration  of  fluoride 
in  water  leaving  those  nodes.   The  computer  model  predicts  a  mixing  zone 
running  just  east  of  Main  Street  along  Garfield,  Beach,  and  Forest  Park, 
while  the  field  tests  indicate  that  fluoride  reached  Main  Street  at  Green 
causing  a  large  mixing  zone  bordered  by  Salem  and  Mishawam  in  northwest 
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Table  21 .   Validation  Test  Fluoride  Concentrations  (ppm) 


Location  (PM  Test)    Sample  No. 


Well  E 

#2  Arlington 

Booster  Pump  (High) 

Elm  and  West 

Rag  Rock  Tank 

#14  Bedford 

Mishawum  and  Garden 

Main  and  Clinton 

Cedar  and  Washington 

#14  Cedar 

New  Boston  and 

New  Industrial 
#204  Salem 
Booster  Pump  (Low) 
#66  Garfield 
#55  Bedford 
#369  Washington 

(near  Olympia) 
#44  Beach  (near  Middle) 
#51  School 

(near  New  Boston) 
#130  Montvale 

(near  Bow) 
Washington  and 

Montvale 
#89  Mishawum  (near 

Forest  Park) 
#576  Main  (near  Kilby) 
Main  and  Middlesex  Canal 
#82  Beach  (near  Salem) 
#87  Olympia  (near 

Mishawum) 
#2  Merrimac 

(near  School) 
#15  Wyman  (near  Main) 
#460  Wildwood  (near 

Olympia) 
Well  F 
Forest  Park 
#21  Clinton 
#12  Alfred  (near  Main) 
#110  Pearl  (near 

Webster  Ct.) 
#9  Fowle  (near  Main) 
#58  Water  (near 

Woburn  Pkwy. ) 
#22  Bow 

#10  Montvale  (near  Main) 
Garfield  and  Green 


1A 
2A 
3A 
4A 
5A 
6A 
7A 
8A 
9A 
10A 

1B 
2B 
3B 
4B 
5B 

6B 
7B 

8B 

9B 

10B 


5C 


Lawrence 

UMASS 

Electrode 

Electrode 

SPADNS 

0.058 

0.06 

0.0 

0.064 

0.07 

0.24 

0.059 

0.07 

0.24 

0.468 

0.48 

0.56/0.40 

0.138 

0.14 

0.30/0.63 

0.088 

0.08 

0.00/0.24 

0.936 

0.98 

1  .34 

0.230 

0.23 

0.62 

0.976 

1  .05 

1  .34 

1  .010 

1 .05/0.98 

1  .26 

0.160 

0.16 

0.24/0.31 

0.966 

1  .00 

0.94 

0.080 

0.07 

0.16 

0.994 

1  .05 

1  .34 

0.148 

0.15 

0.31 

0.610 

0.66 

0.79 

0.551 

0.53 

0.70/0.50 

0.810 


1  .010 


1  .42 


0.170 


0.80 


0.97 


1  .50 


0.16 


1  .2 


1  .2 


>1  .4 


1C 

0.854 

0.82 

1  .2 

2C 

0.195 

0.18 

0.4/0.1 

3C 

0.162 

0.17/0.17 

0.56 

4C 

0.792 

0.76 

1  .40 

0.31 


6C 

0.554 

0.58 

1  .18 

7C 

0.164 

0.16 

1 .26/0.18 

1D 

0.966 

0.94 

1  .34 

2D 

0.076 

0.06 

0.47/0.0 

3D 

0.380 

0.38 

0.63 

4D 

0.200 

0.21 

0.40/0.18 

5D 

0.216 

0.22 

0.56 

6D 

0.194 

0.19 

0.48/0.26 

3 

0.107 

0.10 

K15 

0.066 

0.07 

X18 

0.978 

1  .10 

C24 

0.314 

0.35 

45 

0.970 

1.13 
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Table  21  (continued) 

Montvale  and  Campbell  X46  0.910  0.92 
N.  Warren  and  Harrison  A55  0.066  0.08 
#51  Fowle  (near 

Garfield)  B55         0.954        1  .04 

Montvale  and  Green  X57  0.974  1.02 
Central  and  Utica       66  1.00         1.04 

#69  Wood  67  1 .02         1 .06 

#15  Auburn  (near 

Montvale)  68  1.00         1.10 

#108  Arlington  (near 

Hudson) 
Pleasant  and  Ellis 
#55  Salem  (near  Wade) 
#310  Main  (near  High) 
#7  Green  (near  Main) 
Salem  and  Main 
#80  Lexington  (near 

Well  E) 
#5  Harrison  (near  Winn) 
Main  and  Pleasant 
#117  Green  (near 

James  Terr. ) 
Montvale  and  Washington 
Washington  and  Pine 


D69 

0.085 

0.07 

70 

0.060 

0.07 

73 

0.434 

0.51 

X75 

0.316 

0.37 

89 

0.962 

1  .06 

97 

0.220 

0.25 

SI  31 

0.060 

0.06 

183 

0.077 

0.10 

252 

0.080 

0.08 

A748 

0.978 

1.03 

766 

1  .03 

1  .05 

800 

0.954 
(AM  Tests) 

1  .10 

S15 

0.926 

1  .08 

19 

0.464 

0.44 

S48 

0.994 

1  .00 

59 

0.488 

0.50 

D225 

0.056 

0.07 

Montvale  and  Campbell 
#10  Montvale  (near  Main) 
Montvale  and  Washington 
#438  Main  (near  Everett) 
N.Warren  and  Harrison 
#107  Winn  (near 
Middlesex)  242         0.460        0.55 
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Whispering 
Hill  Tank 


Horn  Pond 
Reservoir 


MDC 


Zion  Hill  Tank 


Figure  17 
Fluoride  (ppm) 
Test:  1-5  PM  6/4/85 
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Whispering 
Hill  Tank 


Horn  Pond 
Reservoir 


Zion  Hill  Tank 


Figure  18 

Fluoride  (ppm) 

Test:  6-7  AM  6/4/85 
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Whispering 
Hill  Tank 


Horn  Pond 
Reservoir 


MDC 


Zion  Hill  Tank 


Figure  19 
Fluoride  (ppm) 

Re-prediction 
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Woburn.   Thus  the  measured  MDC  water  zone  is  slightly  larger  than  that 

predicted  by  the  computer  model. 

V.   Water  Distribution  from  Wells  G  and  H 

The  computer  model  of  the  Woburn  1983-85  water  system,  now  calibrated 
and  tested,  was  changed  to  reflect  the  historic  development  of  the  city.   In 
order  to  go  back  to  the  1 970- ' 79  period,  pipes  48,  83,  84,  and  85  were 
deleted;  nodes  37,  59,  and  60  were  deleted;  pipes  52  and  80  were  changed  to 
six  inches;  and  wells  G  and  H  were  included  at  node  45.   These  changes 
amount  to  a  small  overall  change  in  the  pipe  and  demand  distributions,  but  a 
major  change  in  the  east  Woburn  water  source,  from  the  MDC  connection  to 
wells  G  and  H.   Going  back  to  the  1964-1969  period  required  substantial 
additional  changes  in  the  computer  model.   Pipes  70,  78,  79,  and  81  were 
deleted;  nodes  56  and  57  were  deleted;  pipes  73,  74,  77,80,  82,  and  86  were 
changed  to  six  inches;  and  pipe  80  was  lengthened  to  3500  feet  and  connected 
between  nodes  54  and  58.   These  changes  have  a  large  effect  on  the  northeast 
section  of  Woburn  but  not  on  the  overall  water  supply  system.   In 
particular,  the  user  demand  changes  resulting  from  this  node  removal  are 
quite  small.   Further  user  demand  changes,  needed  to  match  the  user  demands 
of  the  representative  water  distributions  from  wells  G  and  H,  are  made  later 
with  the  scaling  factor. 

The  1970-1979  computer  model  was  compared  to  a  set  of  ISO  test  data 
taken  at  various  times  during  the  years  1970,  1972,  1974,  and  1975  in  order 
to  check  the  validity  of  this  modified  computer  model  and  re-calibrate  the 
model  if  that  proved  necessary.   Table  22  presents  a  comparison  of  the  1970- 
1979  computer  model  and  that  ISO  data.   The  model  used  a  scaling  factor  of 
1.4  as  was  previously  used  with  the  1983-1985  computer  model,  since  the  flow 
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Table  22.      1970-1979   Model   Recalibration 

Node  Data 

No.  H    (ft) 

3  301.7 

3  299.1 

3  304.3 

3  300.6 

6-7  248.5 

21  243.3 

30  254.6 

30-34  250.0 

32  243.5 

42  232.1 

44  240.0 

44  240.0 

50  249.8 

50  238.5 

55  244.9 

55  226.5 
55-61  239.1 
55-61  227.5 

56  226.5 

57  234.1 
57  246.6 
57  246.6 
57-59  240.0 
61  243.1 
68  232.0 
68  211.5 

N  =   26  Avg.   =   3.6   ft    (1 .6   psi) 

S  =   9.5    ft    (4.1    psi) 


Model 

AH 

H    (ft) 

(ft) 

311.0 

9.3 

311  .0 

11.9 

311.0 

6.7 

311  .0 

10.4 

232.3 

-16.2 

244.7 

1.4 

242.1 

-12.5 

241.1 

-8.9 

241  .7 

-1.8 

243.0 

10.9 

243.3 

3.3 

243.3 

3.4 

242.6 

-7.6 

242.6 

4.1 

242.1 

-2.8 

242.1 

15.6 

242.0 

2.9 

242.0 

14.5 

242.0 

15.5 

242.0 

8.5 

242.6 

-4.0 

242.6 

-4.0 

242.6 

2.6 

242.0 

-1  .1 

237.1 

5.1 

237.1 

25.6 
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conditions  of  the  ISO  data  are  not  known.   The  results  show  a  mean  error  of 
3.6  feet  in  total  head,  with  a  root-mean-square  error  of  9.5  feet.   These 
represent  1.6  and  4.1  psi  respectively  and  are  quite  good  since  the  ISO  data 
represents  the  average  of  flow  conditions  over  many  days  in  four  years. 

The  1964-1969  computer  model  was  compared  to  two  sets  of  ISO 
measurements  taken  on  November  3  and  4,  1965  and  on  September  11,  1969. 
Table  23  presents  a  comparison  between  this  data  and  its  computer 
prediction,  again  using  a  peaking  factor  of  1.4.   These  results  indicate  a 
mean  error  of  3.1  feet  in  total  head  but  with  a  root-mean-square  error  of 
16.3  feet,  corresponding  to  1.3  and  7.1  psi  respectively.   The  7.1  psi  error 
is  quite  high,  but  since  the  operating  conditions  and  user  demands  during 
the  tests  were  unknown,  corrective  action  lacks  direction.   Perhaps  the  Rag 
Rock  tank  being  off  line  during  the  September  11,  1965  measurements  is  the 
cause,  but  there  is  just  not  enough  data  to  know. 

The  pumping  history  of  the  wells  and  the  daily  demand  cycle  of  Waldorf 
and  Cleary  (6)  are  not  used  to  describe  the  water  distributions  from  wells  G 
and  H.  Their  eight  months  do  not  represent  the  pumping  history  of  wells  G 
and  H  correctly,  and  their  demand  cycle  is  also  wrong. 

First,  the  peaking  factors  used  in  this  previous  study  do  not  account 
correctly  for  the  daily  operating  cycle  of  reservoirs  filling  during  the 
night  and  emptying  during  the  day.   Their  daily  cycle  was  based  on  the  cycle 
of  the  MDC  inflow  rather  than  user  demand.   Herein  minimum,  average,  and 
maximum  user  demands  are  used  for  each  typical  flow  distribution.   These 
demands  are  based  on  the  constant  pumping  rates  of  the  wells  and  the  inflow 
and  outflow  of  the  city's  reservoirs.   The  Shaker  Glen  booster  pump  is  set 
"on"  for  the  maximum  and  average  parts  of  the  daily  demand  cycle  since  the 
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Table  23. 

1964-1' 

}69   Model 

Recalibrati 

Data 

Data 

Model 

AH 

H(ft) 

Year 

H(ft) 

(ft) 

237.4 

1965 

232.4 

-5.0 

216.9 

1965 

240.1 

23.2 

207.6 

1965 

246.2 

38.6 

246.7 

1965 

247.2 

0.5 

246.4 

1965 

245.2 

-1  .2 

248.1 

1965 

240.7 

-7.4 

240.8 

1965 

245.3 

4.5 

278.8 

1969 

244.2 

-34.6 

272.5 

1969 

243.5 

-29.0 

240.8 

1969 

237.2 

-3.6 

234.1 

1965 

241  .4 

7.3 

251  .7 

1969 

241  .4 

-10.3 

235.5 

1969 

235.5 

0.0 

252.2 

1965 

237.4 

-14.0 

258.0 

1969 

237.4 

-20.6 

235.9 

1965 

235.6 

-0.3 

232.6 

1969 

235.6 

3.0 

256.3 

1969 

234.0 

-22.3 

238.5 

1965 

244.0 

5.5 

244.2 

1965 

245.0 

0.8 

240.2 

1965 

241  .2 

1  .0 

Node 
No. 

6 

7 

9 

18 
22 
34 
43 
44 
51 

53 
54 

54 
61 

63 
63 
66 
66 

69 
23-24 

20-21 
39 

N  =  21  Avg.  =  -3.1  (-1 .3  psi) 

S  =  16.3  (  7.1  psi) 
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Zion  Hill  reservoir's  water  level  charts  show  that  the  pump  was  off  only 
during  the  midnight  to  dawn  period  when  user  demand  was  low.   The  water 
level  charts  for  two  of  Waldorf  and  Cleary's  eight  typical  months  are 
available  and  are  presented  in  Appendix  VI.   They  generally  show  a  maximum 
demand  roughly  20$  higher  than  average  during  six  morning  hours,  average 
demand  during  12  afternoon  and  evening  hours,  and  minimum  demand  roughly  2055 
below  average  during  six  late  night  hours. 

Second,  the  monthly  average  pumping  rates  used  by  Waldorf  and  Cleary 
give  improper  water  distributions  for  wells  G  and  H  since  neither  pump  rates 
nor  flow  patterns  average  in  a  meaningful  way.   Wells  are  on  or  off.   Water 
from  wells  G  and  H  does  or  does  not  reach  a  certain  street.   Thus  three  of 
Waldorf  and  Cleary's  typical  months  should  have  been  divided  into  two  parts 
each;  December,  197^  needs  one  part  when  well  G  is  on  alone  and  the  other 
when  G  and  H  are  pumping  together,  while  November  1967  and  July  1971  need 
one  part  when  well  G  is  on  and  another  when  it  is  off.   Fortunately  the 
on/off  status  of  only  a  few  pumps  are  important.   Due  to  the  Rag  Rock  and 
Horn  Pond  reservoirs,  the  pumping  status  of  wells  A„,  B,  Cp,  D  and  F  has 

little  effect  on  the  flow  distribution  of  Woburn's  water.   Only  the  status 
of  wells  E,  G,  and  H  and  the  Shaker  Glen  booster  pumps  is  important. 

Finally,  the  monthly  average  pumping  rates  used  by  Waldorf  and  Cleary 
were  not  correct  for  those  months  when  wells  G  and/or  H  did  not  pump 
throughout  each  day  but  only  pumped  during  parts  of  some  days.   The  Woburn 
pumping  records  indicate  both  daily  and  monthly  volumes  of  water  pumped  by 
each  well,  the  days  when  pumping  occurred,  and  the  hours  for  each  day  and 
month.   Waldorf  and  Cleary  used  pump  flow  rates  based  on  the  monthly  volumes 
and  number  of  days,  assuming  full  2H   hour  days.   This  method  gives  flow 
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rates  equal  to  the  actual  if  the  days  were  full  days  but  less  than  actual 
when  the  pumps  were  on  for  partial  days. 

Since  the  pumping  history  of  Waldorf  and  Cleary  has  the  deficiencies 
indicated  above,  a  completely  new  history  was  developed.   The  daily  pumping 
records  for  all  days  when  either  wells  G  or  H  were  on  were  examined.   It  was 
noted  that  the  volume  of  water  pumped  was  not  metered  until  April,  1972  but 
was  estimated  by  the  water  system  operators.   Pumping  rates  based  on  the 
volume  of  water  and  hours  pumped  were  calculated.   The  number  of  days  each 
month  that  wells  E,  G,  and  H  pumped  in  their  various  combinations  were 
calculated  based  on  the  hours  pumped  each  day  that  month.  The  monthly 
average  flow  rate  of  wells  A  ,  B,  C„,  D,  and  F  was  calculated  for  each 

month,  as  well  as  the  percentage  of  the  total  monthly  supply  that  came  from 
wells  G  and/or  H.   All  these  figures  are  presented  for  the  1964-1969  and 
1970-1979  periods  in  Tables  24  and  25  respectively. 

The  flow  patterns  for  the  water  from  wells  G  and  H  were  divided  into 
two  groups,  one  with  well  E  on  and  the  other  with  it  off,  since  the  status 
of  well  E  caused  significant  differences  between  those  flow  patterns.   Then 
the  effect  of  the  percentage  of  the  total  water  supply  provided  by  wells  G 
and/or  H  was  examined  and  also  found  to  be  significant.   The  number  of  days 
when  each  percentage  from  wells  G  and  H  occurred  was  calculated  for  each 
historical  period,  1964-1969  and  1970-1979,  and  for  each  well  E  status,  on 
or  off.   Figures  20  through  23  present  the  results  of  those  calculations. 
The  distributions  shown  in  these  figures  are  quite  varied  and  so  the  maximum 
percentage,  the  minimum  percentage,  and  the  median  percentage  were  chosen  as 
representative  of  the  various  percentages  that  occurred  during  the  pumping 
history  of  wells  G  and  H. 
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Table  24.  Wells  G  and  H  Pumping  History:  1964-1969 

Year  &   Mo.  Avg.   Days  Daily     %   from  %   from      Days  Days     Daily 

Month    MGD:  A     Pumped  MGD       Wells  G,  Wells  G,    Well  E  Well  E   MGD 

Wells  G,  Well  G,    H,G  &  H  H,G,  &  H    on  off      Well  E 

B,C  ,D,F   H,G  &  H  H,G  &  H    (E  on)  (E  off)      (G/H  on)  (G/H  on) 


1964 


Oct 

2.80 

30.8 

Nov 

2.63 

23.5 

Dec 

2.54 

26.8 

1965 

Dec 

2.15 

29.9 

1966 

Jan 

2.15 

31  .0 

Feb 

2.33 

28.0 

Mar 

2. .60 

31.0 

Apr 

3.42 

30.0 

May 

3.24 

29.3 

June 

3.80 

30.0 

July 

4.38 

31  .0 

Aug 

4.24 

•31 .0 

Sept 

2.97 

30.0 

Oct 

2.68 

31 .0 

Nov 

3.06 

17.0 

Dec 

2.62 

25.5 

1967 

Jan 

2.68 

25.5 

Feb 

2.72 

28.0 

Mar 

2.80 

12.0 

July 

4.43 

3.0 

Aug 

3.65 

8.0 
20.0 

Sept 

3.18 

4.0 

3.0 

23.0 

Oct 

2.74 

31 .0 

Nov 

2.91 

25.0 

Dec 

4.00 

4.3 

1968 

Apr 

3.83 

19.0 

May 

3.59 

31.0 

June 

3.43 

30.0 

July 

4.37 

31 .0 

Aug 

4.53 

31 .0 

Sept 

4.02 

29.9 

Oct 

3.64 

31.0 

Nov 

3.45 

24.8 

0.72  G  16  20  13.2  17.6  0.98 

0.72  G  17  21  9.1  14.4  1.01 

0.72  G  17  22  11.4  15.4  1 .03 

1.21  G  28  36  3.4  26.5  0.90 


36 
33 
29 

23 
24 

21 
19 
19 
25 
27 
25 
28 


27 
27 
26 
12 
14 
32 

27 

14 

35 
39 
37 


23 

24 

25 

21 

20 

22 

24 

25 


1.25  G 

29 

1.15  G 

26 

1  .05  G 

- 

1  .00  G 

- 

1  .00  G 

- 

1.00  G 

18 

1  .00  G 

15 

1  .00  G 

16 

1  .00  G 

20 

1  .00  G 

22 

1  .00  G 

20 

1.00  G 

22 

1  .00  G 

21 

1  .00  G 

21 

1  .00  G 

21 

0.58  H 

- 

0.58  H 

- 

1  .73  G&H 

29 

1  .15  G 

23 

0.58  H 

13 

1  .7  3  G&H 

31 

1 .73  G&H 

34 

1  .73  G&H 

33 

1.73  G&H 

27 

1  .15  G 

20 

1  .15  G 

21 

1  .15  G 

22 

1.15  G 

19 

1  .15  G 

18 

1.15  G 

20 

1  .15  G 

21 

1.15  G 

22 

3.7 

27.3 

0.93 

2.3 

25.7 

0.86 

- 

31 .0 

- 

- 

30.0 

- 

- 

29.3 

- 

2.7 

27.3 

0.85 

18.9 

12.1 

1  .16 

10.4 

20.6 

1  .03 

15.2 

14.8 

0.95 

14.1 

16.9 

0.97 

8.3 

8.7 

0.93 

9.3 

16.2 

0.98 

9.8 

15.7 

1  .05 

15.0 

13.0 

1.03 

4.0 

8.0 

0.97 

- 

3-.0 

- 

- 

8.0 

- 

5.0 

15.0 

0.67 

3.7 

0.3 

0.67 

2.6 

0.4 

0.67 

12.7 

10.3 

0.67 

22.0 

9.0 

0.66 

23.0 

2.0 

0.64 

4.3 

— ■ 

0.65 

15.0 

4.0 

0.70 

27.0 

4.0 

0.66 

29.9 

0.1 

0.68 

27.1 

3.9 

0.69 

27.4 

3.6 

0.68 

24.6 

5.3 

0.6U 

22.9 

8.1 

0.63 

21.0 

3.8 

0.63 
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Table   24-    (continued) 


Dec 


3.37 


25.5 


1  .15    G 


22 


25 


19.5 


6.0 


0.68 


1969 


Jan 

3.66 

25.2 

1  .15  G 

21 

24 

14.8 

10.4 

0.69 

Feb 

3.48 

28.0 

1.15  G 

22 

25 

15.0 

13-0 

0.69 

Mar 

3.^2 

31  .0 

1  .15  G 

22 

25 

18.6 

12.4 

0.69 

Apr 

3.45 

30.0 

1.15  G 

22 

25 

16.9 

13.1 

0.69 

May 

3.65 

31  .0 

1  .15  G 

21 

24 

19.9 

11.1 

0.71 

June 

4.81 

30.0 

1.15  G 

17 

19 

22.5 

7.5 

0.70 

July 

4.56 

31 .0 

1  .15  G 

18 

20 

24.3 

6.7 

0.74 

Aug 

4.28 

31 .0 

1.15  G 

19 

21 

26.9 

4.1 

0.72 

Sept 

3.73 

30.0 

1  .15  G 

21 

24 

21  .1 

8.9 

0.71 

Oct 

4.16 

0.5 

1  .15  G 

19 

- 

0.5 

■- 

0.66 
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Table  25. 

Wells  G  and  H  Pump 

ing  History: 

1970-1979 

Year  & 

Mo.  Avg. 

Days 

Daily 

%    from 

%    from 

Days 

Days 

Daily 

Month 

MGD:  A 

Pumped 

MGD 

Wells  G, 

Wells  G, 

Well  E 

Well  E 

MGD 

Wells  G, 

Well  G, 

H,G  &  H 

H,G,  &  H 

on 

off 

Well  E 

B,C2,D,F 

H,G  &  H 

H,G  &  H 

(E  on) 

(E  off) 

(G/H  on) 

(G/H  on) 

1970 

Apr 

3.91 

15.0 

1.15  G 

20 

23 

10.6 

4.4 

0.70 

May 

3.78 

21.8 

1.15  G 

20 

23 

17.3 

4.5 

0.72 

June 

4.72 

30.0 

1  .15  G 

18 

20 

25.6 

4.4 

0.71 

July 

4.47 

30.5 

1.15  G 

18 

20 

22.0 

8.5 

0.68 

Aug 

4.29 

31.0 

1  .15  G 

19 

21 

24.4 

6.6 

0.66 

Sept 

3.51 

30.0 

1.15  G 

22 

25 

17.0 

13.0 

0.67 

Oct 

3.90 

15.5 

1  .15  G 

20 

23 

12.4 

3.1 

0.67 

Nov 

3.55 

9.8 

1.15  G 

21 

24 

7.3 

2.5 

0.65 

Dec 

3.47 

10.3 

1  .15  G 

22 

25 

8.0 

2.3  ' 

0.66 

1971 

j 

May 

4.01 

3.1 

1  .15  G 

20 

- 

3.1 

- 

0.68 

July 

4.93 

22.8 

1.15  G 

17 

19 

20.0 

2.8 

0.71 

Aug 

3.73 

31.0 

1  .15  G 

21 

24 

29.3 

1  .7 

0.70 

Sept 

3.60 

27.9 

1.15  G 

21 

24 

15.9 

12.0 

0.68 

Oct 

3.04 

25.5 

1  .15  G 

23 

27 

19.6 

5.9 

0.71 

Nov 

3.02 

16.5 

1.15  G 

24 

28 

16.0 

0.5 

0.67 

Dec 

2.69 

25.2 

1  .15  G 

25 

— 

25.2 

- 

0.67 

1972 

Jan 

2.97 

25.8 

1  .15  G 

24 

28 

25.7 

0.1 

0.67 

Feb 

2.81 

29.0 

1.15  G 

25 

29 

16.7 

12.3 

0.68 

Mar 

2.95 

28.3 

1  .15  G 

24 

28 

17.3 

1  1.0 

0.69 

Apr 

3.30 

16.2 

1  .06  G 

21 

- 

16.2 

— 

0.70 

1974 

• 

* 

Aug 

4.51 

26.8 

1.10  G 

17 

- 

26.8 

- 

0.69 

3.3 

1 .74  G&H 

25 

28 

3.1 

0.2 

0.69 

Sept 

2.98 

2.3 

1  .02  G 

22 

- 

2.3 

- 

0.69 

0.3 

0.58  H 

14 

- 

0.3 

- 

0.68 

27.4 

1 .60  G&H 

30 

- 

27.4 

- 

0.68 

Oct 

2.81 

0.9 

1  .01  G 

22 

- 

0.9 

- 

0.73 

0.3 

0.58  H 

14 

- 

0.3 

- 

0.73 

29.8 

1.59  G&H 

31 

- 

29.8 

- 

0.73 

Nov 

2.78 

10.7 

1  .03  G 

23 

- 

10.7 

- 

0.63 

0.3 

0.58  H 

15 

- 

0.3 

- 

0.63 

18.6 

1.61  G&H 

32 

- 

18.6 

- 

0.63 

Dec 

3.22 

8.4 

1  .04  G 

21 

- 

8.4 

- 

0.64 

4.1 

1  .62  G&H 

30 

- 

4.1 

- 

0.64 

1975 

May 

4.31 

14.1 

1  .07  G 

18 

20 

14.3 

0.1 

0.64 

June 

4.18 

27.7 

1  .04  G 

18 

- 

27.7 

- 

0.63 

88 


Table    25    (continued) 


2.3 

1 .62   G&H 

25 

- 

2.3 

- 

0.63 

July 

4.16 

1.3 

1 .01    G 

17 

- 

1.3 

- 

0.63 

29.7 

1 .59   G&H 

25 

- 

29.7 

- 

0.63 

Aug 

3.79 

13.1 

1  .01    G 

19 

- 

13.1 

- 

0.64 

17.9 

1 .59    G&H 

26 

- 

17.9 

- 

0.64 

Sept 

3.19 

20.7 

1  .01    G 

21 

- 

20.7 

- 

0.64 

0.3 

0.58   H 

13 

- 

0.3 

- 

0.64 

9.0 

1.59    G&H 

29 

- 

9.0 

- 

0.64 

Oct 

3.64 

19.4 

1  .01    G 

19 

- 

19.4 

- 

0.65 

Nov 

4.00 

12.5 

1.04    G 

18 

- 

12.5 

- 

0.65 

Dec 

3.42 

29.8 

1  .02    G 

20 

— 

29.8 

- 

0.66 

1976 

Jan 

3.23 

31  .0 

1  .05    G 

21 

- 

31.0 

- 

0.64 

Feb 

3.69 

28.0 

1  .06   G 

20 

- 

28.0 

- 

0.65 

Mar 

4.01 

17.2 

1  .02    G 

18 

20 

17.1 

0.1 

0.62 

June 

5.24 

0.1 

0.97   G 

14 

- 

0.1 

- 

0.64 

0.8 

0.53   H 

8 

- 

0.8 

- 

0.64 

20.2 

1.50    G&H 

20 

- 

20.2 

- 

0.64 

July 

4/66 

31  .0 

1  .  60    G&H 

23 

- 

31.0 

- 

0.61 

Aug 

4.53 

9.6 

0.57   H 

10 

- 

9.6 

- 

0.58 

21  .4 

1 .58    G&H 

24 

26 

20.8 

0.6 

0.58 

Sept 

4.53 

30.0 

0.57   H 

10 

11 

29.9 

0.1 

0.58 

Oct 

4.38 

•30.0 

0.57   H 

10 

— 

30.0 

— 

0.57 

1977 

Jan 

4.31 

9.5 

1  .10    G 

18 

- 

9.5 

- 

0.58 

Feb 

4.11 

27.2 

1  .06    G 

18 

- 

27.2 

- 

0.56 

Mar 

4.05 

27.8 

1  .06    G 

19 

- 

27.8 

- 

0.56 

Apr 

4.02 

28.9 

1  .05    G 

19 

21 

13-6 

15.3 

0.56 

May 

3.60 

19.6 

1.03    G 

19 

22 

16.6 

3.0 

0.92 

11.4 

1.67    G&H 

27 

- 

11.4 

- 

0.92 

June 

3.86 

27.7 

1  .00    G 

17 

21 

27.0 

0.7 

0.92 

2.3 

1.65    G&H 

26 

- 

2.3 

- 

0.92 

July 

4.. 11 

30.8 

1  .00    G 

17 

20 

30.5 

0.3 

0.91 

Aug 

3.85 

30.9 

0.98   G 

17 

20 

29.8 

1.1 

0.89 

Sept 

3.79 

30.0 

0.96   G 

17 

20 

29.8 

0.2 

0.85 

Oct 

3.64 

31.0 

0.95    G 

17 

20 

29.0 

2.0 

0.86 

Nov 

3.61 

30.0 

0.95   G 

17 

21 

22.5 

7.5 

0.87 

Dec 

4.33 

31.0 

0.97   G 

— 

18 

- 

31.0 

— 

1978 

Jan 

4.42 

31.0 

0.98   G 

- 

18 

- 

31.0 

- 

Feb 

4.44 

28.0 

0.98   G 

- 

18 

- 

28.0 

- 

Mar 

4.55 

5.6 

0.98   G 

- 

18 

- 

5.6 

- 

18.2 

0.69  H 

- 

13 

- 

18.2 

- 

7.2 

1 .58   G&H 

- 

26 

- 

7.2 

- 

Apr 

4.26 

6.0 

0.62  H 

- 

13 

- 

6.0 

- 

24.0 

1 .62   G&H 

- 

28 

- 

24.0 

- 

May 

4.31 

31  .0 

1.66   G&H 

24 

28 

0.4 

30.6 

0.99 

June 

4.30 

29.8 

1  .62   G&H 

24 

- 

9.8 

- 

0.88 

July 

4.45 

31.0 

1 .63   G&H 

23 

- 

1.0 

- 

0.90 

Aug 

4.32 

31  .0 

1 .64   G&H 

24 

28 

8.9 

22.2 

0.81 

89 


Table    25    (continued) 


Sept 

4.43 

0.6 

1  .02 

G 

- 

23 

- 

0.6 

- 

29.4 

1  .62 

G&H 

- 

27 

- 

29.4 

- 

Oct 

3.96 

0.3 

1  .02 

G 

- 

20 

- 

0.3 

- 

30.7 

1  .61 

G&H 

- 

29 

- 

30.7 

- 

Nov 

3.98 

0.1 

1  .02 

G 

- 

20 

- 

0.1 

- 

29.9 

1  .62 

G&H 

- 

29 

- 

29.9 

- 

Dec 

4.31 

31.0 

1.62 

G&H 

24 

27 

20.3 

10.7 

0.77 

1979 

Jan 

4.12 

0.1 

0.59 

H 

- 

13 

- 

0.1 

- 

30.9 

1.61 

G&H 

25 

28 

9.6 

21  .2 

0.80 

Feb 

4.50 

28.0 

1  .65 

G&H 

24 

27 

1.5 

26.5 

0.79 

Mar 

4.06 

31.0 

1.58 

G&H 

24 

- 

1.0 

- 

0.82 

Apr 

3.73 

30.0 

1.55 

G&H 

26 

29 

9.0 

1.0 

0.79 

May 

4.00 

21  .7 

1.55 

G&H 

25 

- 

1  .7 

- 

0.66 
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Figure  20 
1964-1969 
Well  E  on 
Total=  619.0  Days 
Avg.%  =  21.3 
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Figure  21 
1964-1969 
Well  E  off 
Total  =  543.8  Days 
Avg.%  =  26.4 
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Figure  22 

1970  - 1979 

Well  E  on 

Total  =  1359.3  Days 

Avg.%  =20.9 
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Figure  23 

1970-1979 

Well  E   off 

Total  =  473.8    Days 

Avg.%  =  24.1 
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The  boundaries  between  the  part  of  Woburn  served  by  wells  G  and  H  and 
the  rest  of  Woburn  were  calculated  using  the  computer  model  for  each 
particular  historical  flow  distribution.   The  1964-1969  boundaries  are  given 
by  Figures  24  through  29  while  the  1970-1979  boundaries  are  presented  in 
Figures  30  through  35.   In  each  figure  the  solid  line  represents  the 
boundary  determined  by  the  average  user  demand  (peaking  factor  =  1.0)  of 
that  particular  day,  while  the  dotted  and  dashed  lines  represent  the 
boundaries  determined  by  the  maximum  (peaking  factor  =  1.2)  and  minimum 
(peaking  factor  =  0.8)  user  demands  respectively  of  that  same  day.   Thus 
each  figure  shows  the  daily  cycle  of  movement  of  the  boundary  of  the  water 
from  wells  G  and  H.   The  computer  output  corresponding  to  these  figures  is 
in  Appendix  VII. 

Each  of  these  boundaries  divides  Woburn  into  three  zones:   Zone  A  with 
water  from  the  wells  around  Horn  Pond,  Zone  B  with  water  from  wells  G  and  H, 
and  Zone  C  with  a  mixture  of  water  from  both  sources.   This  mixture  occurs 
at  pipe  junctions  along  the  boundaries  of  Zones  A  and  B.   In  many  cases 
unmixed  water  flows  into  a  boundary  node  via  two  pipes  and  then  mixed  water 
flows  out  of  the  node  through  a  third  pipe.   In  this  case  the  boundary 
between  Zones  A  and  B  is  along  the  pipe  with  the  mixed  water  and  the  users 
along  that  boundary  pipe  form  a  very  thin  zone  C.   Further,  in  some  cases, 
usually  in  northwest  Woburn,  the  boundary  pipe  carries  mixed  water  to  a 
larger  Zone  C  which  is  so  indicated  in  the  water  distribution  figures.   At 
each  mixture  node  the  fraction  of  the  mixture  that  comes  from  wells  G  and  H 
can  be  calculated  from  the  results  given  in  Appendix  VII. 

Both  Figures  24  and  30  show  that  mixed  water  reached  the  Rag  Rock 
storage  tank.   However  this  tank  empties  only  into  the  transmission  line 
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Figure  24 
1964-1969 
Maximum  (34%) 
Well  E  on 
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Figure  25 
1964  -  1969 
Median    (21%) 
Well  E  on 
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Figure  26 
1964  - 1969 
Minimum  (13%) 
Well  E  on 
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Figure  27 
1964-1969 
Maximum  (39%) 
Well  E  off 
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Figure  28 
1964  - 1969 
Median  (25%) 
Well  E   off 
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Figure  29 
1964-1969 
Minimum  (12%) 
Well  E  off 
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Figure  30 
1970  - 1979 
Maximum  (32%) 
Well  E  on 
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Figure  31 
1970-1979 
Median  (20%) 
Well  E  on 


103 


Whispering 
Hill  Tank 


Horn  Pond 
Reservoir 


Zion  Hill  Tank 


Figure  32 
1970-1979 

Minimum  (8%) 
Well  E  on 
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Figure  33 
1970-1979 
Maximum  (29%) 
Well  E  off 
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Figure  34 
1970-1979 
Median  (27%) 
Well  E   off 
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Figure  35 
1970-1979 
Minimum  (11%) 
Well  E  off 


10  7 


leading  northward  to  Main  and  Eaton  Streets.   Due  to  the  lack  of  an 
interconnection  on  Bedford  Road  this  flow  is  separated  from  the  east-west 
flows  of  mixed  water  along  Bedford  road. 

The  user  area  at  Salem  and  Washington  Streets,  shown  in  Figure  5  as 
node  44,  is  the  only  user  demand  area  which  received  undiluted  well  G  or  H 
water  whenever  those  wells  were  pumping.   Figures  32  and  35  indicate  these 
minimum-G/H-distribution  flow  boundaries.   Residents,  commerce,  and  industry 
in  this  area  received  more  water  from  wells  G  and  H  than  other  Woburn  users; 
their  total  contact  time  of  2995.9  days  includes  all  the  days  indicated  in 
Figures  20  through  23. 

A  much  larger  section  of  Woburn  never  received  water  from  wells  G  or  H. 
The  user  zones  nearest  Horn  Pond,  that  is  nodes  9,  11,  17,  18,  19,  20,  21, 
22,  30,  and  31,  were  always  provided  with  Horn  Pond  water.   Again  Figure  5 
indicates  the  location  and  extent  of  these  user  areas. 

The  rest  of  Woburn  received  some  water  from  wells  G  and  H,  but  less 
than  the  maximum  amount  due  to  the  2995.9  days  of  pumping.   The  amount 
depends  on  the  four  sets  of  flow  conditions  described  by  Figures  20,  21,  22, 
and  23.   The  median  boundary  lines  of  Figure  25  separate  the  user  zones  with 
more  or  less  than  half  the  619.0  days  of  contact  with  water  from  G  or  H 
under  the  conditions  listed  in  Figure  20.   Figures  28,  31,  and  34  describe 
the  median  user  zones  for  the  conditions  of  Figures  21,  22,  and  23, 
respectively.   The  user  areas  inside  all  four  of  these  Zone  B  boundaries 
received  G  or  H  watc,  undiluted  by  Horn  Pond  water,  on  more  than  1498  days. 
Thus  the  user  areas  indicated  in  Figure  5  for  nodes  39,  41,  42,  44,  49,  50, 
51  ,  and  58  received  G  or  H  water  on  more  than  half  the  days  those  wells  were 
pumping. 
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The  calculation  of  a  lower  limit  on  the  number  of  days  of  G  or  H 
contact  at  a  particular  node  is  straightforward.   For  example,  consider  the 
user  area  at  Salem  and  Main  Streets,  node  23.   During  the  1964-1969  period, 
with  well  E  on,  node  23  got  Horn  Pond  water  for  both  the  minimum  and  median 
fractions  of  G  and  H  water,  but  was  on  the  boundary  for  the  maximum 
fraction.   Figure  20  indicates  that  on  at  least  22  days  a  mixture  of  Horn 
Pond  and  wells  G  and  H  water  reached  user  area  23.   The  composition  of  this 
mixture  can  be  calculated  from  the  inflow  data  to  node  23  that  is  given  in 
Appendix  VII.   In  the  same  period,  with  well  E  off,  node  23  was  supplied  by 
Horn  Pond  for  the  minimum  fraction  of  G  and  H  water,  and  received  a  boundary 
mixture  on  at  least  the  nine  days  corresponding  to  the  maximum  fraction. 
However  this  time  the  daily  peaking  factor  affected  the  source  of  the  water 
reaching  node  23  for  the  median  fraction.   Node  23  was  supplied  by  Horn  Pond 
during  the  18  hours  of  average  or  greater  user  demand  but  received  a 
boundary  mixture  during  the  six  hours  of  less  than  average  demand.   One 
quarter  of  the  271 .9  days  corresponding  to  flows  with  the  median  fraction  or 
more  G  and  H  water  is  68.0  days,  indicating  at  least  that  many  nights  of 
boundary  mixture.   Thus  the  total  number  of  boundary  mixture  days  at  node  23 
during  1964-1969  is  at  least  22+  9  x  (3/4)  +  68.0  =96,8  days.   Similarly, 
the  1970-1979  period  has  at  least  (3/4)  x  18.6  +  (1/4)  x  73.9  =  32.4  days  of 
boundary  mixture  at  node  23.   Thus  node  23  is  never  within  zone  B,  receiving 
undiluted  water  from  wells  G  and  H.   However  the  user  area  of  node  23  does 
receive  boundary  mixture  water  on  at  least  129.2  days. 

A  similar  calculation  indicates  an  upper  limit  at  node  23  of  1218.4 
days. 
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VI.   Conclusion 

The  purpose  of  this  study  is  to  evaluate  and  complement  the  water 
distribution  study  of  Waldorf  and  Cleary  (6).   Therefore  the  similarities 
between  their  report  and  this  one  are  now  reviewed,  and  the  differences  in 
the  data,  method  of  analysis,  and  results  are  emphasized. 

The  pipe  networks  were  basically  the  same,  with  roughly  100  pipes  and 
70  pipe  junctions.   Roughly  10  percent  of  the  networks  were  different  with  a 
few  important  differences.   Four  of  these  are  the  closed  valve  on  Burlington 
near  Lexington,  the  lack  of  connection  between  the  16  inch  and  ten  inch 
pipes  on  Bedford,  the  six  inch  versus  12  inch  pipes  on  Salem  and  the  pipe 
roughness.   The  user  demand  distributions  were  also  similar,  but  Waldorf  and 
Cleary  did  not  address  the  uncertainty  due  to  the  substantial  pipe  leakage 
in  the  system,  whereas  this  study  uses  a  demand  scaling  factor  of  1.4  to 
convert  the  identified  (metered  plus  residential  estimate)  user  demand  of 
1984  to  the  total  demand  indicated  by  the  calibration's  pumping  records. 
They  also  used  a  peaking  factor  nearer  1 .0  in  describing  the  daily  cycle  of 
user  demand  as  well  as  an  incorrect  monthly  pumping  rate  history.   Finally, 
Waldorf  and  Cleary  did  not  use  the  ISO  data. 

Both  studies  used  pipe  network  analysis  computer  programs  to  transform 
the  data  into  flow  distribution  patterns.   The  methods  differed  in  that  this 
study  used  the  ISO  reports  and  the  related  water  system  data  for  November  2 
and  3,  1983  to  calibrate  both  the  demand  scaling  factor  and  pipe  roughness. 
The  field  tests  used  to  verify  the  calibrated  model  were  not  done  by  Waldorf 
and  Cleary  since  they  restricted  their  attention  to  the  period  when  wells  G 
and  H  were  in  operation.   They  did  not  compare  their  model  with  the  1964- 
1979  ISO  data  either. 
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Even  with  all  these  differences  the  results  are  roughly  the  same.   Both 
studies  indicate  a  boundary  between  the  G  and  H  water  and  the  Horn  Pond 
water  that  runs  from  Montvale  and  Green  northwest  to  Main  and  Mountain  for 
the  1964-1969  period.   The  present  study  indicates  mixed  water  in  northwest 
Woburn  during  the  1970-1979  period  while  Waldorf  and  Cleary  show  Horn  Pond 
water  there.   However  the  present  study  indicates  that  the  variation  in  the 
location  of  this  boundary  is  much  greater  than  that  reported  by  Waldorf  and 
Cleary.   The  wider  variation  is  most  evident  along  Bedford  Road  where  this 
report  indicates  that  a  mixture  of  water  from  the  two  sources  sometimes  even 
reached  Shaker  Glen. 

The  similar  results  of  the  two  studies  indicate  that  the  error  level  of 
the  boundary  location  is  lower  than  would  be  expected  from  the  lack  of  data. 
The  principal  areas  of  ignorance  were  the  user  demand  and  leakage 
distributions,  the  pipe  roughness,  and  the  effects  of  aging  on  the  well  pump 
characteristics.   Fortunately  the  fluoride  distribution  part  of  the 
validation  tests,  shown  in  Figure  17,  can  be  compared  to  the  computer  models 
results  (Figure  19)  for  that  test  in  order  to  get  a  good  indication  of  the 
error  level  in  the  fluoride  boundary  location.   These  figures  show  that  the 
error  level  of  the  boundary  between  Horn  Pond  and  MDC  waters  is  between  zero 
and  one  nodes  since  the  model  and  the  measurement  show  either  the  same 
boundary  nodes  or  adjacent  boundary  nodes.   For  examples,  Salem  and  Beach 
streets  define  an  important  boundary  node  in  both  model  and  field  test  while 
Green  and  Garfield  is  the  model's  boundary  node  located  next  to  the  measured 
boundary  node  at  Green  and  Main  Streets.   This  example  of  boundary  level 
error  of  between  zero  and  one  nodes  is  a  good  estimate  of  the  overall 
boundary  location  error  in  Figures  24-35. 
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